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ABSTRACT: The tensor network/geometry correspondence is a proposed discrete 
version of the holographic duality. We show how important features in the AdS/CFT 
dictionary, such as the bulk operator reconstruction via the HKLL relation and the 
map between bulk isometry and boundary global symmetry, can emerge naturally 
from the tensor network construction. Furthermore, we propose that the tensor 
network living on the Bruhat-Tits tree gives a concrete realization of the recently 
proposed p-adic AdS/CFT (a holographic duality based on the p-adic number field 
Q,); in particular, the wavefunction of the tensor network defines the ground state 
of the boundary p-adic CFT. 
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1 Introduction 


The principle of holographic duality states that quantum gravity in a spacetime M is 
equivalent to a quantum field theory on the boundary 0M [1, 2]. Both conceptually 
and mathematically, it provides one of the most powerful tools to understand non- 
perturbative quantum gravity (via its dual quantum field theory description). The 
most prominent and tractable holographic duality arises when the spacetime is anti- 
de Sitter (AdS), where we have the AdS/CFT correspondence [3]: quantum gravity 
in an asymptotically AdS spacetime is equivalent to a conformal field theory (CFT) 
living on the boundary of the AdS spacetime. 


While the duality was first engineered from string theory, it is believed that the 
holographic principle is independent from string theory. Particularly when consid- 
ering generalized examples beyond those that follow directly from string theory, the 
AdS/CFT dictionary was guessed based on global symmetries and various physical 
and mathematical intuitions. It is by no means obvious how and why, the dic- 
tionary emerges in this correspondence. This is an extremely difficult question to 
answer, because the AdS/CFT correspondence, in most of the parameter regime, is 
a strong/weak duality. 


The study of entanglement in AdS/CFT has provided new food for thought. In 
addition to other aspects of the AdS/CFT dictionary, there is now an extra entry 
— namely the Ryu-Takayanagi (RT) formula — that calls for explanation. The 
RT formula is an important characteristic of semi-classical gravity [4]. It raised the 
questions of the meaning of semi-classicality and locality. 


These pieces of the AdS/CFT dictionary serve as much as hints as they are 
puzzles to our understanding of quantum gravity. The RT formula inspired the 
observation of a suggestive similarity between the tensor network construction of 
wavefunction and the AdS/CFT correspondence [5]. The observation stimulated a 
surge of investigations [6-26] and it appears that the similarities between these two 
different realms of physics might be much more than similarities — the tensor network 
is perhaps capturing some essential physics of the AdS/CFT correspondence. 


In particular, it has led to recent proposals of the holographic code based on 
perfect tensors on a network embedded in negatively-curved space [19]. It was also 
found that perfect tensors emerge very naturally in spaces of tensors with large bond 
dimensions [20]. It was found that such tensor networks naturally recover the RT 
formula and recover some features expected to follow from causality that mimic an 
error correcting code. It was also observed that close to a background of perfect 
tensors, other features of the AdS/CFT dictionary, such as structures that behave 
like Witten diagrams in the computation of correlation functions, also emerge [21]. 


Given all the appealing features of the tensor network, there are some urgent 
conceptual questions. A tensor network lives on a discrete space. To which extent 
can it capture the AdS/CFT correspondence, which, strictly speaking, has only been 
defined for continuous spacetime? Does there exist a discrete version of AdS/CFT 
correspondence which tensor network can capture fully? 


With the goal of answering these questions eventually, in this paper we focus on 
two important aspects of holographic correspondence — bulk operator reconstruction 
and the interplay between bulk isometry and boundary global symmetry — and ask 
whether, and if yes how, they manifest themselves in the tensor network construction. 


In AdS/CFT, a local field in the bulk can be reconstructed from local gauge 
invariant operators on the boundary, convoluted using the HKLL (Hamilton-Kabat- 
Lifschytz-Lowe) kernel (called “smearing function”) and integrated over the bulk 
causal patch [27]. In the tensor network, (using the operator pushing) we find quite 
generically that a discrete version of the HKLL relation [27] naturally appears. More- 
over, the analogues of 1/N corrections to the HKLL relation [28] also emerge nat- 
urally in the tensor network. Let us note that the (linear) HKLL relation, admits 
an interpretation of a wavelet transform. We were inspired by the discussion in [24], 
where the exact holographic mapping was based on the Haar wavelet, to discuss the 
HKLL relation both for the continuous and also the p-adic AdS/CFT in this light. 


In the traditional, continuous, AdS/CFT correspondence there is a mapping be- 
tween bulk isometry and boundary global symmetry. The discrete space on which 
the tensor network lives also has certain isometries. We then demonstrate how ge- 
ometric isometry of the tensor network recovers global space (time) symmetries of 
the wavefunctions. And we show how an interesting subtlety in the cutoff surface 


allows dynamical information regarding scaling dimensions of primary operators to 
be retrieved from the rigid geometric setup. 


In both questions above, we encounter a problem of how to best harness the 
isometry of the tensor network. In the bulk operator reconstruction, the wavelet 
should be decomposable as eigenmodes of the isometry of the network. For the 
boundary symmetry recovery, to make more quantitative use of symmetries, it is 
necessary to have control of explicit wavefunctions living on the tensor network that 
carry well-defined quantum numbers under the geometric isometries. 


However, thus far, our understanding of isometries of tensor network has been 
rather primitive.’ The main problem is that the discretization that the tensor net- 
work realizes are usually chosen to be the dual graph following from some regular 
tiling of the AdS space. Each such tiling preserves a discrete subgroup of the isom- 


etry of AdS (i.e. SL(2,R) in a 2d bulk) generically called a Coxeter group. To our 
knowledge, the representation theory is not very much discussed in the mathematics 
literature, nor are the characteristic functions with well-defined eigenvalues that live 
in the corresponding lattice. 


To gain more symmetry and also quantitative control of graph wavefunctions so 
that we could examine in detail bulk reconstruction from the tensor network, we turn 
to a simpler, but perhaps less physical, construction, namely the p-adic AdS/CFT. 
Inspired by recent discussions [26, 29, 30], we consider tensor network living on a 
tree — more specifically the Bruhat-Tits tree, which is a geometrical presentation 
of the p-adic expansion of a p-adic number. The tree preserves the full conformal 


group SL(2, Qp) — only with the real field R replaced by a different field completion of 


rational field Q. This is much larger than the discrete subgroup of SL(2, R) preserved 


by any regular tessellation of the AdS space. 


The tensor network based on the Bruhat-Tits (BT) tree then gives a concrete 
realization of the p-adic AdS/CFT correspondence proposed in [26, 29]. These papers 
generalize the AdS/CFT dictionary to the situation where the boundary theory is 
taken to live in a space-time that is based on the field Q, (which is continuous), 
and the discrete BT tree plays the role of the bulk AdS space.” Besides being an 
interesting holographic duality that is based on Q, instead of R, p-adic AdS/CFT 
might also become relevant to the analysis of the continuous version via the adelic 


construction. 

In the current paper, we see that features of the p-adic AdS/CFT, such as trans- 
mutation of bulk isometry into boundary global symmetry and bulk operator recon- 
struction, naturally emerge from the tensor network living on the BT tree, which 


‘See however [21] that initiated the study on symmetries in the tensor network. 

?The discussion is mainly based on a two dimensional bulk and a one dimensional boundary. 
The proposed higher-dimensional generalization involves finite algebraic extensions of Q,, for more 
detail see [26, 29]. 


gives us a concrete realization of the p-adic AdS/CFT. At least as far as a bulk 
scalar theory is concerned, it emerges in the tensor network and that its relation 
with the boundary theory follows the usual rules of the AdS/CFT correspondence 
as proposed in [26, 29]. The correspondence between tensor network and a bulk 
geometry certainly deserves further detailed study and test. 


The paper is organized as follows. In section 2 we briefly review the basics 
of the tensor network construction of wavefunctions. In section 3 we explain how 
analogues of the HKLL relation and its 1/N corrections follows from the tensor 
network construction. In section 4 we discuss the relationship between geometric 
symmetries of the tensor network and the boundary wavefunction. Then in section 
5 we review some details of the p-adic numbers, and discuss the tensor network 
embedding in the BT tree. Various details are relegated to the appendices. 


2 ‘Tensor network 


2.1 A lightening review of tensor networks meeting the AdS/CFT 


Solving for the exact ground state of a quantum many-body system analytically is 
in general a very difficult problem, because of the gigantic dimension of the Hilbert 
space. Consider a many body Hamiltonian H. The goal is to find its ground state 
|). Generically, a wavefunction is a rank-N tensor, where N = the number of the 
constituents in the many-body system. i.e. 


W) = XO finin lin +i), (2.1) 


iein 


for 1 < ip < Dk, where D, is the dimension of the Hilbert space at site k. Here, we 
have implicitly assumed that the many body problem admits a Hilbert space that 
can be factorized into direct products of spaces. Determining the wavefunction would 
therefore involve determining [ [F D, number of unknowns, which scales exponen- 
tially with NV. It has long been known in the literature that the entanglement entropy 
of low lying eigenstates, and particularly the ground states of local Hamiltonian sat- 
isfies an area law [31] — the entanglement entropy of a region in configuration space 
scales as the area of the boundary of the region in the large region size limit. In 
the case of 1+1 dimensional critical systems, there is a logarithmic violation of the 
area law. These observations inspired the introduction of tensor network as a (nu- 
merically) efficient way to approximate a wavefunction. Tensor networks are clever 
re-organization of the wavefunction f;j,...;,,, where it is expressed as the contraction 
of other tensors. 

facin = > Taa apase TT (2.2) 
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In the above equation, [; 9... denotes different choices of tensors, each with multiple 
indices. Some of these indices 7; to iy are the physical indices that remain un- 
contracted, whereas un denote auxiliary indices that are contracted between tensors. 
These contractions are most conveniently expressed graphically where vertices denote 
tensors, and legs contracted indices. 


(a) 


(b) 


Figure 1. These are popular tensor networks describing (ground state) wavefunctions. 
(a) MPSs network specialized in 1+1d gapped systems. (b) MERA network specialized in 
1+1d gapless systems. (c) PEPs networks often chosen for 2+1d systems. Picture courtesy 
of [32]. 


Figure 1 gives a few popular architectures of these tensor networks. The graph 
is chosen depending on the entanglement structure. For example, MPS (please refer 
to figure 1) are frequently used for gapped systems, whereas MERA is a popular pick 
for gapless systems in one dimension with build-in logarithmic violation of the area 
law [33]. The design of MERA is worth some discussion. As one can see in figure 
1 (b), the MERA network is made up of layers of tensors starting from the bottom 
layer containing the physical indices. Each layer is meant to be a linear map that 
maps the system to a coarse grained one. Therefore an extra dimension develops as 
coarse graining progresses up the graph. 

At first sight the introduction of these auxiliary tensors would increase the com- 
plexity of the problem. However, imagine that each tensor has P elements to be 
determined, and the graph contains altogether L nodes. The determination of the 
wavefunction would then require fixing L x P elements. If L changes slowly with the 
number of site N — in the case of MPS L = N — the complexity of the problem is 
greatly reduced. The simplification allows the use of feasible numerical methods to 
determine the values of the constituent tensors by minimizing the expectation value 


of the energy (2)|H|w). More details about tensor networks can be found for example 
in [32, 34]. 

Tensor network architects appeared to live in a very different world from string 
theorists until very recently. It took a prescient observation of Swingle [5] to con- 
nect the tensor network and the AdS/CFT correspondence. It was pointed out that 
the MERA type tensor network, as alluded to above reproduces some basic features 
of holography — the entanglement entropy of a subregion is bounded from above 
by the length of some geodesic cutting through the network, which bears a strong 
resemblance to the Ryu-Takayanagi (RT) [35, 36] formula. As already mentioned 
above, the idea of the MERA tensor network encodes the process of coarse graining, 
such that an extra dimension develops in the graph, and this is not unlike our in- 
terpretation of the extra radial direction in the AdS bulk which plays the role of an 
RG scale. This led to a natural conjecture: that the tensor network encapsulates the 
underlying physics of the AdS/CFT correspondence. The suitable choice of a graph 
that defines the tensor network is in fact a discretization of the (AdS) geometry. 
Since the graph is directly encoding the wavefunction of a many-body system, the 
correspondence between geometry and state is exact. 


There are various progresses that render the correspondence between tensor net- 
works and some interacting theories perhaps containing gravity in AdS space more 
plausible and concrete. For example, it is observed that in addition to the RT for- 
mula, the bulk field reconstruction using boundary operators in the AdS/CFT dic- 
tionary is consistent with a network that carries error correcting properties [13, 14]. 
These observations provide intuitions into possible algebraic properties of tensors 
that would actually resemble semi-classical AdS geometry and a concrete proposal 
is made in [19]. There, networks defined on a negatively curved graph with con- 
stituent tensors given by “Perfect Tensor” (Ty, ,9,--a.,) are found to saturate the RT 
formula exactly. Perfect tensors are even-rank tensors that satisfy the following two 
properties [19], 


e For any partition of the 2n indices into two groups J = {a1,---a,} and a = 
{Qn41,°**@2n} with equal number of indices, Tja is a unitary map: 


Tial op = Ong (23) 

e For a partition in which the number of indices in J is less than that in a, Tja 
is a norm-preserving projection operator from a to J: 

Taly = zy DË, (2.4) 


where D is the bond dimension i.e. the number of values taken by each index 
of the tensors, and s counts the difference between the number of indices be- 
tween the collection J and a. Roughly speaking, the bond dimension D plays 


the role of the rank N of the boundary gauge group in traditional AdS/CFT 
correspondences. 


These properties are crucial in establishing the RT formula. Interested readers are 
referred to the original work [19] for the details of the derivation. Related proposals 
are developed in [20, 37]. The analysis of the random tensors in [20, 38] establishes 
that perfect tensors emerge when the network is composed of typical tensors of large 
bond dimension, supporting the notion that a semi-classical limit emerges in some 
kind of large N limit. 


There are also alternative proposals that discuss the connection between MERA 
and the kinematic space and also works that consider continuous MERA and the 
AdS spaces [8-12]. 


2.2 Open questions in tensor network 


There are many important open questions to be addressed if the correspondence 
between tensor networks and the AdS space is more than an analogy. We enumerate 
in the following a few most glaring ones. 


e Bulk operator reconstruction. 


The perfect holographic code construction was partially inspired by the pattern 
of bulk operator reconstruction in the AdS/CFT dictionary according to the 
HKLL formula [27], in which a bulk operator can be recovered using only a 
subregion in the boundary [13, 14]. The perfect holographic code however 
is found to produce no connected correlation function between local operators 
[21]. It turns out, as we demonstrate in detail in the appendix, that this implies 
that a bulk operator can be reconstructed only as macroscopic products of 
boundary operators. This is in sharp contrast with the HKLL formula where 
in the large N limit, the leading contribution is linear in the boundary operator. 


e Interplay between global symmetries of the wavefunction and isometries of the 


tensor network. 


As it is well known in the AdS/CFT correspondence, global space-time sym- 
metries of the boundary CFT manifests themselves as isometries in the bulk. 
This was one of the most important entry in the AdS/CFT dictionary. Any 
conserved currents must also imply the existence of corresponding gauge fields 
in the bulk. This has to be quantitatively reproduced and harnessed if the 
tensor network captures the correct physics. Some discussion that the tensor 
network does reproduce such features already appeared in [39]. 


e Time dependent evolution of the wavefunctions and dynamics in the bulk. 


While a wavefunction can be evolved using any Hamiltonian that acts on the 
boundary legs, it is not clear whether we can interpret that evolution as some 
local dynamics of bulk degrees of freedom. 


In this paper, we will take some further steps in clarifying the first two questions. 
As we already briefly mentioned in the introduction, both the first and second item 
have a natural interpretation in the tensor network, even at a quantitative level, 
although we are only scratching the surface of these issues. Question three remains 
largely open and we hope to deal with it in a future publication. 


Before we end this section, we note that it is also discussed extensively in the 
literature that the perfect code produces the wrong entanglement spectrum. There 
are various proposals that attempt to resolve this issue based on introducing a center 
in the operator algebra [41]. We took the viewpoint that perhaps this is not as serious 
an obstacle as it initially appeared. The reason is that the MERA type tensor 
networks have been in use in explicit numerical applications in which gapless ground 
states are simulated, and entanglement entropies obtained to very good accuracy. 
In any actual system, the bulk tensors are never given by perfect tensors, and such 
properties are only expected to emerge approximately when the bond dimension is 
large. Therefore, the problem is perhaps not so serious an issue if we were working 
with an explicit Hamiltonian and when we solve honestly for the bulk tensors say by 
numerical methods. We hope to revisit this issue in the future. 


3 Bulk operator reconstruction 


The holographic code is generically an isometry map from boundary degrees of free- 
dom to bulk degrees of freedom. Therefore operators acting in the bulk can generi- 
cally be mapped to operators in the boundary. Recall for example the pentagon code 
in [19], for each bulk operator acting on the bulk vertex, it can be replaced by opera- 
tors acting locally on the legs of that vertex, making use of the isometric property of 
the tensor. E.g. X — XX XXX in the pentagon code. In turn, operators acting on 
the legs of these tensors can be further pushed down to operators acting on the other 
legs of the neighboring tensors. This process can be reiterated until all the operators 
are acting on the boundary dangling legs. This lies at the heart of the bulk operator 
reconstruction using operators within the causal wedge, a property first noticed in 
the AdS/CFT which takes the form of the HKLL operator reconstruction [27]. 
There is however an obvious mismatch between the perfect pentagon code and 
the HKLL relation. In the pentagon code, we see that a bulk operator is equivalent 
to a product of local operators, rather than a sum of local operators acting on 
individual legs at the boundary. In fact we can readily show that for perfect code, it 
is not possible that an operator be represented as an operator acting on fewer than 
half of the legs. We will defer the proof to the appendix. This is indeed an alternative 


way of reiterating the observation in [21] that perfect codes do not support connected 
correlation functions between local operators. 

Suppose therefore that we are dealing with a tensor network with tensors devi- 
ating from perfection. Therefore generically the operator acting on the bulk index 
can now be mapped to action on (a subset of) its legs, and similarly action on a 
given leg is equivalent to the action on the other legs of the same tensor. This is 
precisely the meaning of operator pushing as it is defined in [19], except that we 
are now working with non-perfect tensors. This condition suffices to preserve the 
error correcting property implied by causality in the AdS space. To solve for these 
relations explicitly, we assume that we divide the indices of the tensor T at a node 
into two groups. If there are an even number of indices say 2L, we separate them 
into two groups each of L indices. Then we assume that this matrix, as a linear 
map from L indices to the other L indices is a non-degenerate matrix that admits 
an inverse T~'. Compared to the perfect code, we have thus relaxed the condition 
to make T unitary in every way of dividing the indices, but only non-degenerate. If 
the total number of indices in a tensor is odd = 2L+1, then we place L tensors on 
one side and L + 1 on the other side which we are going to push the operator to. 
In this case, inverse corresponds to contracting the L + 1 indices and recovering the 
identity operator on the other L indices. 3 With these rules the result for operator 
pushing takes the following general form (We are illustrating an example where each 
tensor is a 3-leg tensor) : 


1, (On eo Lutte = Soa niu) (Oj Jaa! & Opp’ F OF (kw, iw) Odal & (OF ow 
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(3.1) 
where we consider the linear map between operator acting on the leg k,, belonging to 
a tensor at node w and operators acting on other legs iw, Jw of the same tensor. The 
left hand side of equation (3.1) is illustrated in figure 2. In general there is no reason 


why a‘ should be identical to a. But to ensure that the pushing of an operator 


to different legs give the same result, we would like to put a¥ = ã&¥ as a constraint. 
The J’s enumerate the operator basis acting on the link. Here we consider the subset 
of O* which can be equivalently understood as an operator acting on one of the 
legs of the tensor. This is possible because in a given node at which tensors T/ ,, ... 
reside, the bulk space J and edge space li, l2,--- have the same dimension and are 
related by a unitary map. i.e. D|Z] = []/ D[l,], where Dfa] denote the dimension of 


the indices a, and L is the total number of legs the tensor possess. 


3In this case, there is an extra subtlety that the inverse is not unique, and we can construct a 
projection operator P to project the rhs of (3.1) to the L leg subspace. The projection operator for 
the 3-leg tensor for example, takes the form Pab, ar = (D>, Tip Liew): The coefficients in (3.1) 
are solved in the presence of the projection operator conjugating the rhs of (3.1). 


Figure 2. This illustrates the left hand side of equation (3.1). 


In an isometric tensor network (in which all directions are the same), there are 
constraints over the allowed choice of the tensors to ensure that these coefficients 
a have at least the same eigenvalues independently of which direction we choose 
to push our operator. This is analogous to the constraints based on the transfer 
matrices that were defined in [21]. These constraints and the relationship with the 
transfer matrices will also be relegated to the appendix. 

Repeating this pushing process all the way to the boundary, we would generically 
obtain such an analogous expression as above, which however now involves a sum of 
operators at the boundary 


OF keu) —> > Af g,(w)O5, Eg ` AF audv(W)O5,O dy ap ee (3.2) 


where Qu, Ov now denote the position of physical sites located at the boundary of 
the tensor network. 


3.1 Bulk operator reconstruction — Linear contribution 


Let us focus on the linear contribution of a local boundary operator to a given bulk 
operator. We will also first focus on the linear term in (3.1) and delay the discussion 
of the role played by the non-linear terms in (3.1). Their role will become transparent 
very soon. In a general graph then, an operator that is acting on a particular leg 
adjacent to a node w can be mapped to an operator on a different leg of the node w 
using (3.1). Subsequently, since the leg is contracted with another node w1, we can 
further choose to push the operator across node w so that it acts on yet another 
leg of wı. This process can be continued until it hits the boundary. The collection 
of a product of a} (kw, iw) as we push the operator along a path that connects the 
starting bulk point w to a boundary point Ov gives a contribution to the coefficient 
Af ay: In a general graph, the complete contribution from the linear term in (3.1) to 
AF o» is thus a sum over products of af (kw,,tw,) along all the paths that connect w 
and Ov : 

AT 5y(W) limsat = ` I] af (ey) Ein); (3.3) 
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where P,, are paths connecting w and O,, and e,, are edges along the path. 

Now focus on the product. Suppose this is a homogeneous network in which 
ay .(v) has no v nor e dependence. Then we can diagonalize this matrix a‘, with 
eigenvalues Às. In this basis, the linear term in the boundary reconstruction of 


eigen-operators O(A) would be given by 


AŤ 5, (w) = OKT ` AS Tapo = prx (Pu lw) = Gog (w, O»), o = 
Pw 


where d is the distance, or number of edges in path P connecting the point w and 
the boundary leg ðv. We have deliberately written A in terms of a new parameter 
A and the valancy p+ 1. We have defined the expression G, (w, ôv). As we will see 
below, in the case of a tree graph this is precisely a graph Green’s function of the 
graph Laplacian (of fixed valency p + 1) describing propagation between w and Ov. 
We will also briefly discuss the case with variable valency. Our construction bears a 
strong resemblance to path-integrals in which we sum over some kind of amplitudes 
across different paths. There are many constraints that have to be imposed on a¥ 
in order for the emerging theory to make sense as a quantum theory in its own right. 
For example, we would expect that the propagator dies off with increasing distances. 
This is only possible if the eigenvalues of a/* is less than 1. There are many more 
issues with propagation around loops, but we will defer a complete classification of 
consistency conditions in a future publication. Let us also note that generally and 
more so for graphs with loops, to give a global definition of the function G,,.(w, v) it 
is necessary to give a global flow analogous to the discussion in [37]. This flow needs 
not be conserved, but rather irrotational. For a given flow, the function G,,.(w, 0») 
so obtained is a green’s function of a directed graph with orientation defined by the 
flow. We will return to this discussion in a future publication. 

To avoid some of these subtleties we will in the following restrict our discussion 
to tree graphs. In that case, the path between any two given nodes on the tree is 
unique. As a result, Az; reduces to the product of a; over the unique path along 
which the operator is pushed down from the bulk node. 


Af a, (w ea. linear — II ak (ev; End (3.5) 
v; EPw 


where Pu now denotes the unique path connecting the bulk node w and the boundary 
dangling leg, and e, denotes the edges emanating from some bulk vertex v that 
belongs to the path P,,. This is illustrated in fig 3. 

In the following, we would like to demonstrate that both (3.3) and its simplified 
version (3.5) for tree graphs are indeed the graph propagator. 


—l1-—- 


Ov 


Ou 


Figure 3. Linear contribution of boundary operators to a bulk operator depends on the 
path that connects them. 


3.1.1 Graph Laplacian and bulk boundary propagator 
A generic graph Laplacian is defined as [43] 
elv) = Sd’) — V(w) o(v), (3.6) 


v~u 


where v’ are all neighbouring vertices of v, and V(v) is the valancy of the given 
node v. In the following however we will first focus on graphs with fixed valancy V 
independent of v. 

3.1.2 (p+ 1) valent graph 


The Laplacian on a V = p+ 1 valent graph can be written as 


=T,- (p+), (3.7) 


where the action of T, is defined as 


(TA) = X oe’), (3.8) 
for all v’ that is the neighbor of v i.e. they are connected by an edge. T, is called 
a Hecke operator. We note that it has the same eigenvalues as the Hecke operators 
defined on the upper half plane with Eisenstein series their eigenfunctions. 

One can show readily that the Green’s function satisfying 


(0, + m?)G(u, V2) = ô(v1, V2) (3.9) 


is given by 
Guy, v2) = N pAder), (3.10) 


=a 


h 
where , 1 r 


m =- oaa = F Tap Pan 


1- p~s 
This is the same relation as discussed in the p-adic BT tree in [29], which is nothing 


other than a tree with valency p+ 1. As promised, the graph Green’s function (3.10) 
is indeed what emerges in (3.5). One can also check that (3.3) would satisfy the same 
equation (3.9) since the Laplacian is a local equation, but would keep the Green’s 
function well defined since d(w, v) is no longer unique in a graph that admits loops. 


3.1.3. Wavelet decomposition vs graph Green’s function 


The discussion above demonstrates that the linear part of the bulk boundary operator 
reconstruction reproduces salient features of the HKLL relation — namely, that it 
naturally is given by some propagator that lives in the emergent geometry, in this 
case the graph on which the tensor network wavefunction is defined. On the other 
hand, in previous works, it has been observed that the holographic mapping from a 
bulk tensor network to a boundary wavefunction is naturally constructed using the 
idea of wavelets — in [24] and more recently in [25], the choice of tensors in the bulk 
is based on the Haar wavelet, where the 3-valent tensor is chosen such that it defines 
an isometry that maps a pair of fermionic operators acting on neighbouring site to 
the symmetric and anti-symmetric basis. 

ast? z Bem + a3,), (3.12) 
where s denotes the layer number and 7 the site number on the layer. These coeffi- 
cients +//2 parametizing the contribution of operators in layer s can be understood 


precisely as the af, in the discussion above. Wavelet transforms thus appear to 
be interpreted as the linear contribution of the generic tensor network. Given this 
observation, it is thus natural also to interpret the HKLL relation as a continuous 
wavelet transform. We will make the statement more precise in the following. 


3.1.4 HKLL bulk reconstrunction as wavelet transform 


The wavelet transform is a technique used frequently for example in signal processing. 
The idea is that depending on the properties of the signal involved, it is frequently 
more convenient to work with basis other than the Fourier modes. In the following, 
we will review briefly some general aspects of wavelet transform, and we will borrow 
the notations in [44], whose discussion makes the parallels with the HKLL relation 
particularly transparent. 

The set of wavelet basis containing d+ 1 parameters is defined as follows. First, 
choose the mother wavelet Ņ(z), which is a (local) function of d parameters 7. The 
set of daughter wavelets is generated from the mother wavelet by translation by @ 
and rescaling by s: 


), (3.13) 
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then form the set of the wavelet basis. 


Given a signal function f(z), the wavelet transform using the wavelet w is then 


defined as i 48 
~ ta = 
WE) = | ate Feu" —)F@), (3.14) 
The inverse transform is given by 
E 1 °° ds 7 = 
f(z)= A 3H [ae Wy (G, s) Pasz) (3.15) 
where ppe 
oo k 2 
o= f" a POE 9 
0 


and (K) is the Fourier transform of o(Z) and k = |k|. For the inverse transform to 
be well-defined, we need the following “admissibility condition” 


Cy < +00 (3.17) 


It is often more convenient to apply the wavelet transform in the Fourier space. 
Consider Fourier transforming (3.14) w.r.t. @. This gives 


> 


W,(R,s) = 9(k,s)f(R) with gf) =~ 


(vt) (=k) (3.18) 
It means that to define the mother wavelet of a wavelet transform, we could equiva- 
lently specify a function g(k,s) which is a function of the momentum and the scale 
factor s. 

In fact, we could have defined the wavelet transform most generically as 


£ Wss) = gAn A, (3.19) 


© where f (1) is the transform of f w.r.t. to some generic orthnormal complete set of 
modes in the space where f(z) is defined: 


f(x) = fa FOxa (x), fe VOX, (2)xa;(@) = AnA), Pxa (x) = Ax(2), 

(3.20) 
for some differential operator D defined in space x and some appropriate measure 
for the eigenvalue space ;. In terms of A; then the wavelet transform is 


W;(a,s) = | deae too, (3.21) 


and moreover 


h= T dds gAn 8)xt, (a) xa (2). (3.22) 


= j4— 


Then all of these can be readily recasted in terms of Fourier modes if we replace 
A, > k, Xa (£) > exp(ikx) and D > 0”. 

Let us pause here and note that the expression in (3.22) described above actually 
appears naturally in many occasions. One such occasion is the heat kernel. The heat 
kernel is defined the solution of the equation 


0,H(#,s|7) +O,H(z,s|7) = 0, (3.23) 


and hence can be written in terms of eigen-modes of the Laplacian 


H(#,s|7) = f Daon (3.24) 


where e~™5 


à plays the role of g(\,s). The heat kernel however does not admit an 
inverse transform as a wavelet, since the expression in (5.52) is not convergent when 
we take W(k) = e™*. In the following, we will explore in further detail the HKLL 


relation in light of the wavelet transform. 


In AdS/CFT correspondence, a normalizable local field in the bulk ¢(%, z) * that 
is dual to the boundary operator O with conformal dimension A, i.e. lim,_,9 (7, z) = 
z7^O(T) can be reconstructed from the boundary operator via 


$(#,2) = / d'y K(#,2|7) OW) 


+E fata! cade.) | din KEDO) | din KE, ZO) 


M | (3.25) 


where K (2, z | y’) is the normalizable solution of the (Lorentzian) bulk equation of mo- 
tion. Usually, to recover the boundary field, the extrapolation dictionary is adopted, 
in which we read off the boundary value of the bulk field $(x,z) to recover the 
boundary operator. Here, given the comparison with wavelet transforms, we hope to 
complete the discussion in an analogous manner by defining an inverse transform. To 
obtain the smearing function K(g,z |y) in configuration space in a compact form, 
an analytic continuation to de-Sitter space is performed in [45] which send spatial 
coordinates x + ix. This gives a smearing function that looks like 


RL] x 


g)? 
here and from now on 7° = 6,72". G(Z,z|z',z’) is the bulk green function, with 
fall off limy G(Z, z|Z’, 2’) ~ 2-4 L (2, z |Z) + 274K (z,z| 2’), where L(Z, z |7") is 
the normalizable solution of the bulk EOM. 


“Here (Z, z) denotes the Poincaré coordinate with the metric ds? = 4 (d2? + njdx'dz’), with z 
being the radial direction and 7 the boundary d-vector. 


d—A 
KE- (ip) e-e- (3.26) 
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The leading term of the bulk reconstruction (3.25) is then very close to a wavelet 
transform with the mother wavelet 


1 d—A 
AL -2 
palT) = (; = =) O(1 — 7°) (3.27) 
where © is the step function with O(x) = 1 for x > 0 and O(x) = 0 for x < 0. Her 
daughter wavelets, defined via (3.13), and the HKLL “smearing function” are related 
by a scaling 

K(#,2|f) = 2° *de2(9) (3.28) 


It is tempting to interpret this result as the fact that the AdS/CFT correspondence 
selects for us a set of wavelet transform: for each operator Oa of conformal dimension 
A, there is a mother wavelet Ya such that the wavelet transform Wo(#, z), which is 
related to the bulk fields ¢(z, z) via 


d 


olz, z) = 2^7 3Wo(Z, z) (3.29) 
becomes a weakly-coupled and semi-classical degree of freedom. The holographic 
direction “z” in the Poincare coordinates is identified with the scaling parameter of 
the wavelet transform. 

Now let’s look at the inverse transform. For the standard inverse wavelet trans- 
form (3.15) to exist, the admissibility condition (3.17) needs to be satisfied. For a 
scalar field, the conformal dimension is A = d/2 + v with v = ,/d?/4+m?. The 
Fourier transform of (3.27) is 


A J (k 
palk) = a (3.30) 
which gives 
Jk? 
Cy = / dk nA) ~ T(0) > œ (3.31) 


i.e. Cy diverges and the inverse wavelet transform (3.15) is no longer valid. However, 
this does not pose a problem for us because, although we use the wavelet trans- 
form to obtain the bulk field, we actually do not need to use the (standard) inverse 
wavelet transform (3.15) to obtain the boundary operator from the bulk field. We 
can compute instead 


f d'z dz ./GAdSa $(2, 2)K (2, z0), (3.32) 


which has an extra dressing factor of z2” relative to the integration measure of the 
(standard) inverse wavelet transform (3.15). (The scaling parameter “s” there is to 
be identified with “z” in the Poincaré coordinate here.) This is a natural choice in 
the AdS geometry. 
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The Fourier transform of the smearing function is [45] 


dj. „d/2 F 
d*k z% TAKZ) ika- 


K(Z, z| =N r) 7 


(3.33) 


time-like 


where k = v —k2 and M the normalization constant. The integral is restricted to 
time-like momenta so that J, remains normalizable. Plugging (3.33) into (3.32), and 


using 
m © dz 1 
p= — Aker = — 3.34 
Gsf Z nts (3.34) 
for time-like momenta k, we are left with 
N? d?k A ao 
— k” Oal’) ie, 3.35 
2v time-like (27)? olz) ° l ) 


Note that the integral (3.34) is the normalization that replaces Cy in the wavelet 
transform (3.15). The normalization condition Ĉ, < +00 in the inverse transform 
(3.32) replaces the “admissibility condition” (3.17) that appears in the wavelet trans- 
form (3.15). The restriction to time-like momenta is due to this normalization con- 
dition, without which the integral is not finite. The restriction however means that 
the reconstruction of the original operator Og is restricted by causality. Namely, 
equation (3.35) can be re-written as 


2 


a d@z'Oa(z')G(z' — Y), (3.36) 
V 
where ji 
G(z# — 7) = kT ett (G-#) 3.37 
( 5) -_ (2n)¢ l ) 


This should be contrasted with the usual wavelet transform where the admissibility 
condition is usually satisfied for all k. 

When we put the tensor network on the p-adic BT tree, these results would again 
naturally emerge there. We will postpone the discussion of the case of the p-adic 
AdS/CFT as another explicit example in a later section. 


3.2 Non-linear contribution 


In the previous subsection, we have focused on the linear contribution in the oper- 
ator pushing map at one site (see equation (3.1)) and discussed its contribution to 
the linear term in the bulk operator reconstruction. Let us generalize this one step 
further, and consider quadratic contributions to the bulk operator. This will eventu- 
ally lead the way towards a complete interpretation of the bulk-boundary operator 
reconstruction in (3.1). 

Given the discussions in the previous sections, a very general picture has emerged. 
The linear contribution in (3.1) defines a free theory on the graph. The non-linear 
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terms in (3.1) then behave as interaction vertices. They contribute to both linear 
terms and non-linear terms in (3.2). 

Let us consider the following. We would like to obtain systematically the val- 
ues of AT auas W) in (3.2) at fixed w, ðu, Ov. Given the discussion on the linear 
contribution, the answer is immediate. 

As we use (3.1) to map operators acting on legs from a node across neighboring 
nodes to their legs, we are allowed to choose to branch out to two operators across 
one nodes i.e. picking up the coefficient OFF iy jy exactly once across some node v. 
Across any other nodes we are only allowed to pick up the contribution of the linear 
term ar, across any other nodes. Recall that when we pick up only af, along all 
paths connecting two nodes u, v, that can be interpreted as the propagator satisfying 
the equation (3.9). As a result, we have 


AP autos (w) = > Gox (w, va 4, 54 Gor (v, OU Gay (v, Ov2), (3.38) 


where the sum is over all possible nodes at which the branching occurs. 

The situation again simplifies when the graph concerned is a tree, in which the 
propagators G simplifies. Moreover, it is not hard to convince oneself that for trees, 
the vertex v at which the splitting occurs is actually unique. It always lies on the 
leg that is crossed by the geodesic that connects the two boundary points v1, Ov». 
This is illustrated in figure 4. 


Ov 


Oavi Oaz 


Figure 4. Non-linear contribution of boundary operators to a bulk operator depends on 
the path that connects them and also an interaction vertex marked blue in the figure. 


We mentioned that these non-linear terms in (3.1) contributes to linear terms 
in (3.2) as well. The relation (3.1) alone is not enough. We need to combine it 
with another set of relations in which multiple operators acting on different legs are 
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combined into a single operator upon pushing across a node : 
OF (ky) O” (lw) = > Big” (Favs bo; My)OM, +++ (3.39) 


With (3.39), we can include extra contributions to A% (w) in which we choose 
a path from node w to some vertex v, at which point it splits to two operators. Then 
these two operators propagate independently but they are joined at another node u 
using (3.39). We call this a 1-loop contribution to A‘: 


AK, itoop = X` Gox lw, v) 0% Go, (0,U)Go,(0,U) BE Gou (u, ðv), (8.40) 
where we have omitted indices indicating edges to avoid clutter. The generalization 
of the above (multi) loop corrections to other non-linear terms at the boundary is 
immediate. We will discuss in the appendix how the splitting into two operators 
across a node is closely related to the fusion matrix that is defined in [21]. This 
discussion recovers precisely the systematics of the 1/N corrections to the HKLL 
relation [28]. For the simpler case of tree tensor networks however there aren’t 
further loop corrections. 


4 Isometry and global symmetry in a tensor network 


In the AdS/CFT correspondence, one very important aspect of the dictionary is that 
global symmetries of the boundary CFT states are recovered as isometries of the bulk 
spacetime. This should emerge in the tensor network construction. This is already 
quite apparent in the earlier works’ discussion on correlation functions of [20, 21], 
where it was shown that the negative curvature of the graph with appropriate geo- 
metric isometries must be responsible for the emergence of correct scaling behaviour 
of correlation functions. 

Relationship between global symmetries of wavefunctions and its manifestation 
in tensor network constructions has been considered in the case of MERA in [39]. It 
was recently reconsidered in the context of projected entangled pair states (PEPS) 
and matrix product states (MPS) [47], where very rich algebraic structures are ob- 
tained. 

In this section, we will discuss the encoding of global symmetries of a wavefunc- 
tion in a tensor network, although some of these results had been known in the tensor 
network literature for sometime. 

We will also demonstrate in two examples how the boundary global symmetry 
manifest as transformations on the bulk tensors — it involves both geometric shifting 
of the tensor sites and the rotation of the tensors. Then we reverse the direction: we 
start from the bulk tensor and ask what geometric symmetry a given tensor network 
can furnish. We focus on the regular tessellation of the Poincare disk. It can furnish 
two different types of symmetry groups: a discrete subgroup of the bulk isometry 
SL(2, R) or the full conformal group SL(2, Qp). 
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4.1 Realizing boundary global symmetry via tensor transformations 


For a state invariant under some symmetry transformation, it means 
gl¥) = |¥). (4.1) 


The state |W) however is constructed from contraction of tensors in a network. This 
means that the description of the state in terms of these tensors involve a large 
amount of redundancy — any contracted leg between two neighboring tensors T7 
and T7 can be rotated by a unitary transformation: 


nets Ties tee T, I E ee (4.2) 


an M 


This transformation U is often also referred to as a gauge symmetry in the standard 
tensor network literature. 


Therefore it is not necessary that individual tensor T stays invariant under the 
transformation g. We only need to ensure that these tensors are invariant up to a 
gauge transformation. In the following, we will discuss the cases of internal symme- 
tries and spacetime symmetries separately. 


4.1.1 Global internal symmetries and bulk gauge symmetries 


For tensors with some physical dangling legs, equation (4.1) implies that 


= T(V) aie: ta Wia i UA (e2) AES (4.3) 


cal 


Jab T(v)o 


“y po 
vey steg? 


where W7 u, (ej) are gauge transformations acting on edge ej, and te, are the indices 
on the link, such that the neighboring tensor contracted to the same edge e; has to 
be transformed by the same conjugate gauge transformation W3 yle): In general 
these W4 (e) need not be the same as g and can certainly admit e dependence. Now 
by starting with the most exterior layer where it is in direct contact with physical 
sites, the tensor replaces the external g transformation by W9 on the virtual legs as 
in (4.3), which in turn requires that the next layer of tensors not in direct contact 
with physical sites behaves as 


Wi it, (C1) Tee, stegs (0) = Ties stare (U) Wit ieg (€2) 0 (4.4) 


in 


The above requirement on the tensor is equivalent to the following: 


Lice’) = ws't Gye il (0) Wi (e2) Dae (4.5) 


: + . 
lej ste, 1 e2? 


where the tensor is exactly invariant under a set of gauge transformations induced 
by the physical global symmetry transformation. As commented in [39, 40], this is 


5We use capital Latin letters to label tensors on individual sites of the network, Greek letters 
for contracted indices, and small Latin letter for boundary legs. 


— 20 — 


the correspondence between global symmetry and gauge symmetries in the tensor 
network. While the tensor network is by construction invariant under local gauge 
transformation (4.2), a global symmetry of the wavefunction dictates that the tensors 
in the interior have to be invariant under a gauge transformation individually. This is 
the analogue of the correspondence between global gauge transformation and global 
symmetry of the boundary theory. Such invariance relation (4.5) defines the invariant 
gauge group (IGG) in an MPS discussed in [47]. Incidentally in MPS the IGG is 
responsible for the projective nature of representations of the transformation at the 
physical sites. Whether the interesting algebraic structure explored in [47] could be 
defined in the MERA type tensor networks is currently work in progress. 


Equations (4.3), (4.4) and (4.5) are based on constructing wavefunctions in- 
variant under global internal symmetries, explaining why the vertex index v stays 
invariant everywhere. 


4.1.2 Spacetime symmetries and graph geometric isometries 

When it comes to spacetime symmetries, the corresponding relation (4.3) should be 
replaced by 

(e1)W? 


a ap 
leg teg 


EN (4.6) 


In a holographic type tensor network, this is very interesting. In addition to 
inducing gauge transformations, the tensor at the top level has moved elsewhere. It 
means that we would also have to replace (4.4) by a similar relation, in which the 


JT (V)o,it, ity =T (9-V)aie ini Wa 


a of 
teiste 


tensor on the R.H.S. of the equation has to move to a different vertex g.v. The 
question is where it should move to. 


The above equation only considers the simple situation in which motion of the 
vertices in the bulk suffices. In general, the symmetry transformation might also 
involve the permutation of the legs of a tensor. This would also require a systematic 
study. Since internal symmetries are relatively well understood compared to space 
time symmetries, we will illustrate our discussion in this subsection with two explicit 
examples. In the following, we will consider the Pentagon code proposed in [19]. We 
will pick suitable graphs to illustrate how both kinds of transformations could come 
into play. 


4.1.3 Examples using the pentagon code 


We now study the example of the pentagon code. We will consider two different 
types of wavefunctions to illustrate the above point. For the 5-qubit pentagon code 
[19] there are two different types of tensors of the form Tege: The I index can take 
two different values and the lower indices takes values of either 0 or 1. The upper 
index is generally referred to as the bulk index. 

Now we consider the following prototype network as our first example: 
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Figure 5. The boundary transformation corresponds to interchange the legs a,b,c to 
a’,b’,c’. Also d,e, f to d’,e’.f’and g,h,i,7 to g',h’,i’,7’. The corresponding bulk trans- 
formation required to keep the wavefunction invariant is to interchange the bulk vertices 
(I,J, K, L) 0 (T', J’, K', L’). 


The corresponding wavefunction is 
pI J T J K K' L L' 
lW) =T cap! def aylarve wT "el fl aly! By'55'eL | 'y6' de! ~ ghij Liwii é! (4.7) 
labedef ghija'b'c d'e f g'i jh 
where |abcdef ghij a'b'c d'e f' g'R'i'j') is the basis ket corresponding to the physical 
degrees of freedom. Let’s look at the following translation operator acting on the 
boundary legs: 


Ti = [[ Ue F +X eX! +Y gY + Z e Z’) (4.8) 


<ij> 


where the superscript i denotes the site index. X,Y, Z are the three Pauli matrices 
and J is the identity matrix. The operator J” exchanges the it site and the jt? 
site. The product runs over all such pairs, with an overall effect of translating every 
bulk index forward. We note that the network chosen above is inspired by the tensor 
network version of the BTZ black hole constructed in [21]. This direction along 
which there is translation symmetry actually corresponds to z > Az and z > Ax 
in Poincare coordinates. The periodic identification ensures that the transformation 
is an exact symmetry of the network even though the network is finite. We will 
therefore call it the scaling transformation albeit it is somewhat misleading here. 
Shortly we will discuss how scaling symmetry is in fact realized in the finite graph. 
Using this operator we will perform a scaling transformation of the boundary legs. 


=D). 


At the boundary, consider the symmetry transformation in which we translate the 
legs a,b,c to a',b,c', d,e, f to d’,e’, f’ and p,q,r,s to p’,q',1r’, s’. This can be achieved 
using a combination of the operators (4.8), giving 


|) =(1 + Xam Doane + Vad alan + Zoazi Zakes) 
I J I! J’ K 7 L L' 
Lite ap! def ya el wp tae fl ary L Bry b6te! Bist Set 1 ghij ed gthiitj é (4.9) 
labedef ghija'b'c d'e f' g/h't'7'), 


where |abc a,b,c a2b2C2 a3b3c3d3) is the basis ket corresponding to the physical de- 
grees of freedom. 

It is not hard to prove that |i) obtained above is precisely equivalent to the 
state obtained by translating the bulk nodes following bulk isometries. Specifically, 
we need to interchange the four bulk sites with unprimed indices with their primed 
partners: (J, J, K, L) © (/', J’, K', L’). This is again achieve using the 7, now acting 
on the bulk indices. To demonstrate that the resultant state is equivalent to (4.9), 
we need to use the stabilizers of the perfect tensors. For example, 


XG Teie = -Tapy de Load po yo (4.10) 
Similarly for the other two Pauli matrices. 

Yrs TS de = -TZ py de Xaa Yeo Xye (4.11) 
and 


Zig ‘ae = —T py de Yaa Zeb Yye (4.12) 


The are more interesting structures involved in the bulk-boundary symmetry corre- 
spondence. To illustrate a different aspect, we consider another example. In this 
example, we find that global transformation on the boundary can be effected in the 
bulk by unwinding of bulk legs in addition to moving all the bulk nodes. Consider 
the following network illustrated in figure 6. 
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Figure 6. The boundary transformation corresponds to interchange the legs a,b,c with 
f,e,d. The corresponding bulk transformation required to keep the wavefunction invariant 
is to interchange the bulk vertex J with K, which are denoted by the green and the red 
nodes respectively, and we have to twist the contracted legs, i.e. we need to interchange 
x,y with z,w. 


The wavefunction corresponding to this tensor network is 
lY) = T” avczy T” zarez T yüjw T jure pig T" zwaes abca'b' ddef d'e f'g'). (4.13) 
We apply the following transformation to the boundary physical legs: 


Ti) =(1 F X anit das F Yam Yaz F Zaia Zad) 
(1 + Xoy Xey + Youn Yey + Zby Zeya) 
(1 F Xe pei fes F Yezi Ee F Zezi Zf) 


I L J M K ted Ith ot 
T EE zabiz T güjul jde fig T nina, abca b'c'defd'e fg Je 


(4.14) 


The above boundary transformation again moves the boundary legs along the peri- 
odic direction. The transformation can again be equivalently achieved by transforma- 
tion in the bulk. In this case however, the new ingredient is that we will have to act 
on the interior legs along with the vertices. To be precise, the bulk transformations 
can be done in two steps. First, we exchange the interior legs: 
W’) =(1 + Ag X zaz + Yogi zas + A ai Zama) 
(Lae Apes F Yun wg F Sane) (4.15) 
T' avery maver: T wel” patil zwaeplabca'b'c'defd'e f'g") : 


Then exactly as we did in the previous example we exchange the bulk vertices, con- 
verting the J, K indices into their primed counterpart. Next we exchange the J and 
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K vertices of |b"). Altogether, we finally have: 


W) = + Xp XK +YrrYk e + Zi ZKx') 


(1 RÄ A ip tes ega + Vict os + VATN Zaza) (4 16) 
(1 al Xy X wn + Yyan Yw, F Zya Zuya) 
T' avery maver T yiiz jwe pig T” waeplabca'b'c'defd'e' f'g’) : 
Now using the operator pushing [19] we can easily show that 
lY) = |"). (4.17) 


These examples illustrate the set of transformations in the bulk that can be con- 
sidered as bulk diffeomorphism. Applied globally they generate boundary global 
symmetry transformations, in a manner consistent with the AdS/CFT dictionary. 


4.2 Bulk isometry group furnished by the graph 


In the previous subsection, we demonstrate that geometric symmetries of the graph 
encodes boundary global spacetime symmetries. Therefore our choice of graphs would 
allow us to select a particular (discrete) subgroup of the conformal symmetry of the 
boundary wavefunction to be preserved explicitly. One systematic way of describing 
AdS isometries in a graph is to embed it in a regular tiling of the hyperbolic space. 
The discrete subgroup preserved by a given tiling is characterized by the Coxeter 
group. A brief review and relevant references can be found in [21]. In [21], the 
symmetries were harnessed to describe orbifolds and a construction of discrete graph 
corresponding to a BTZ black hole. 


This is however only a very crude application of the bulk symmetries. In the usual 
continuous version of the AdS/CFT correspondence, wavefunctions in the AdS space 
can be solved explicitly, and these solutions can be organized into representations of 
the conformal group, allowing us to identify bulk representation of primary operators 
and so on. Our current understanding of the Coxeter group however has not allowed 
us such luxury. We are not aware of a systematic discussion of graph wavefunctions 
defined on regular tilings of the hyperbolic space, or if such solutions exist at all, 
how they organize themselves into representations of the Coxeter group. 

The main reason for this difficulty is that the symmetry preserved by the regular 
tiling is too small: if we insist on the tensor network be geometrically embedded 
in the hyperbolic space and respect the isometry of the substrate, i.e. the isometry 


preserved by the tensor network is a subgroup of SL(2, R), then the remaining isom- 


etry (preserved by the tensor network) is necessarily a discrete subgroup of SL(2, R). 
Most fatally, this discrete subgroup does not contain a discrete subgroup of the scal- 
ing subgroup of SL(2, R). In order to have a better understanding of the quantitative 


interplay of graph symmetries and that of the boundary theory, we need to have a 
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bigger symmetry at our disposal. Interestingly, this can be achieved by relaxing some 
conditions on the construction of the tensor network based on regular tilings. 


We have seen that the isometry of the tensor network based on regular tiling is 


a discrete subgroup of SL(2, R) because we insist on the tensor network be geometri- 


cally embedded in the AdS», respecting its isometry. This is a natural assumption, no 
less because we expect the tensor network to realize a discretize version of AdS/CFT. 
However, there is no a priori reason that the discretization should work in this most 
naive form. Suppose we relax this assumption, i.e. the lattice is first considered as 
an abstract lattice, free from the underlying AdS, then it is possible for the lattice 
to furnish different, possibly bigger, symmetry. 


The second, independent, way to maneuver for more symmetry is to remove 
certain set of edges in the lattice if necessary, thereby give the lattice more freedom 
to have bigger symmetries. (Note that we still choose to place a tensor in every node 
of the lattice given by the regular tiling.) In [21], all nearest neighbors are connected. 
The opposite choice is to draw the least number of lines connecting all sites, i.e. we 
first draw a spanning tree on the lattice,° then only contract those pairs of tensors 
that are connected by the edges of the tree. Therefore a regular tessellation can 
also give rise to a tree-type tensor network.’ This is illustrated in figure 7. The 
tensor network built on the spanning tree would be different from the original tensor 


network on the full lattice, but could retain many essential features.® 


Naively, one might think the spanning tree breaks more symmetry than the 
original lattice. This is true if we demand that they both furnish the discrete part 


of the isometry group SL(2,R) (of the underlying AdS substrate) that is preserved 
by the corresponding tree or lattice. However, since we have also relaxed this latter 
assumption, we end up with an abstract tree, which can a priori have much bigger 
symmetries. 


Since modeling AdS/CFT is our main goal, we still want this symmetry to be 
related to the conformal symmetry in some way. It turns out that for those spanning 
tree with (p + 1) valence, where p is a prime number, it furnishes a representation 
for the full conformal group of SL(2,Q,) where Q, is the field of the p-adic number. 


This is much bigger than a discrete subgroup of SL(2, R). 


6Given a lattice, its spanning tree is a tree that contains all the vertices of the lattice and has 
the minimal number of edges. 

“Even for the same tree, if we make it represent different symmetry groups, the tensor network 
built from it would be different. This ties in with another important ingredient in the bulk boundary 
correspondence: the cutoff surface. Different interpretation of the same tree gives different (natural) 
cut-off surfaces, which in turn determines different boundary dynamics. 

8We leave the study on the relation between the spanning tree tensor network and the lattice 
tensor network, and the relation between tensor networks based on different choices of the spanning 
tree (of the same mother lattice), to future work. 
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Figure 7. The tree graph with valancy =3 embedded in the Poincare disk tiled according 
to the triangular group W (oo, 2,3). 


This consideration is inspired by the recently discussed p-adic AdS/CFT [26, 
29]. The proposed duality is a discrete analogue of the AdS/CFT. In the simplest 
example where the bulk is two dimensional and the boundary one dimension, the bulk 
geometry is given by the Bruhat-Tits (BT) tree, whereas the boundary is purportedly 
a theory that is defined on the field Q, i.e. the p-adic numbers, and which preserves 
the SL(2, Qp) symmetry. (We will review the p-adic AdS/CFT in Sec. 6.) To explain 
the symmetry group SL(2, Qp), we will first review the p-adic numbers and Bruhat- 
Tits tree, a (p + 1) valent tree furnishing the SL(2,Q,) symmetry, as the p-adic 
counterpart of the hyperbolic plane. 


5 The p-adic number field and the Bruhat-Tits tree 


We have seen that for the tensor network to be fully developed the underlying lattice 
(the network) needs to have enough symmetry. Usually the presence of the lattice 
breaks the continuous isometry group down to a discrete subgroup. However, in the 
case of Bruhat-Tits tree, which is a lattice representation of the p-adic number field, 
the full conformal group is preserved. Therefore a tensor network based on Bruhat- 
Tits tree would have much more symmetry than its counterpart living on the regular 
tessellation. In this section we review the p-number field Q, and its Bruhat-Tits tree, 
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to prepare for the discussion on p-adic AdS/CFT and to fix notation. For textbooks 
on p-adic numbers, see [48, 49]. For its applications in string theory or other fields 
of mathematical physics, we recommend [50-56]. 


5.1 From the rational field Q to the field Q, of p-adic numbers 


The field of rational numbers Q can be extended to the field of real numbers R, 


with respect to the Euclidean norm |z|, i.e. R contains all the limit points of Cauchy 


sequences of rational numbers, using the Euclidean norm || as the metric. Another 


way to view the field R is as the set of numbers with all possible decimal expansions: 


N 
ceR = r= Ù a~l” with O<a,<9 (5.1) 


n=— o0 


where it is possible to have infinite expansions in the negative (but not the positive) 
exponent of 10. The Euclidean norm of the real number |x| satisfies a few axioms. 


1. (elo 


2. |z| =0 <= r=0 


æ 


3. |xy| = |x| |y| for all z, y € 


4. |x+y| < |x| + |y| for all x,y € R. (ie. triangle inequality) 


Starting with the field Q, it is possible to extend it in other ways, with respect 
to different norms that obey the above axioms. Given a fixed prime number p, a 
x € Qcan be uniquely expanded in terms of powers of p: 


oO 


t= ` an p” with an € Fp (5.2) 
n=—N 
where F, denotes the residue field consisting of integers 0,1,...,p — 1. Instead of 


measuring the Euclidean length |x|, the p-adic norm |x|, of x measure its congruence 
with respect to the prime p. The definition is as follows. The expansion (5.2) can be 
rewritten as 

Lp = pr?) ` an p” with aọ#0 and v,(x) EZ (5.3) 


n=0 


namely, the leading term in the p-adic expansion of x is p’?‘). The p-adic norm Elg 
is defined as 
Izl =p? with vlr) EZ (5.4) 


and we define |0|, = 0. Namely, the p-adic norm gives the inverse of the highest power 
of p dividing x. (Note that the more divisible x is w.r.t. p, the smaller norm it has.) 
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Take a rational number a = 23 for example, its p-adic norms for p = 2,3,5, 7,11 


are 
24 1 24 1 24 24 


= =i = 5 =a 


bsi (nem 6s 
T= 55 H7 . 


It is then easy to check that the p-adic norm obeys the first three axioms for the 
norm. The fourth axiom (i.e. the triangle inequality) can also be written as 


jz] <1 ==> jl+a2)<C=2 (5.6) 

The p-adic norm corresponds to the stronger inequality with the choice C = 1 instead 
of 2, and gives [57] 

Iz + ylp < max(|z|p, |ylp) (5.7) 

which is called “strong triangle inequality”.° This strong triangle inequality (5.7) 


also implies 
|z + 2p < |2|p (5.8) 


which violates the Archimedes principle |x + «| > |x| — hence the geometry based 
on p-adic norm is called non-Archimedean. 


We thus see that the rational field Q can have infinitely many different norms: 
the Euclidean norm |z| together with the p-adic norms |x|, for each prime p. The 


real field R is only one possible extension of Q, using the Euclidean norm |z|. Now, 


for each prime p, we can have a different extension of Q using the p-adic norms |z|,. 
Given a fixed prime number p, the field Q, consists of all possible formal expansions 
of the form: 


Q, = {Lp = > an p” | an E Fp} (5.9) 
n=—N 


The p-adic norm (5.4), in particular |p”|, = =, ensures that the formal series (5.9) 


0 


1 
pr? 
converges.! 


The p-adic norm is used to define the following subset of Q,, which will be useful 
in the later construction of Bruhat-Tits tree and the discussion on p-adic integration. 
First, the unit sphere in Qp consists of £p with unit norm: 


U, = {£p E Qp | |z| = 1} i.e. gz|p = ao + ap + ap’ +... a9 #0 (5.10) 
The unit ball of Qp is inside U,: 


Zp = {£p E Q,| |x|, < 1} i.e. Zp = ao + ap + agp” + agp... (5.11) 


°Note that the original triangle inequality is trivially satisfied by the p-adic norm: |x + y|p < 
Lely + Luly. 

10Note that in contrast to the decimal expansion for the real number x € R, for a p-adic number 
Lp E Qp, we allow the expansion to be infinite in the direction of the positive exponent of p but 
not along the negative direction, precisely because a higher power of p has a smaller p-adic norm. 
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Note that the unit ball Z, is precisely the ring of p-adic integers. However, unlike 
Z (which is open in R), Z, is both open and closed (“clopen”) in Qp. This will 


be crucial later in the construction of the Bruhat-Tits tree and in viewing it as 
the (holographic) bulk of the Q, boundary. Finally, we denote the set of non-zero 
elements in Q, as Q% = Q,/{0}, which is Q = [[nczp” Up. 


It is quite remarkable that the Euclidean norm |z| and the p-adic norm |z|, 
are the only possible norms to complete the rational field Q (giving R and Q,, 
respectively), as already shown by Ostrowski in 1919 [58]. Conceptually this already 


makes the p-adic fields Q, very interesting: they are the only siblings of R (as possible 
completions of Q). Moreover, Q, is much more than a mathematical curiosity. In 
algebraic number theory questions can be asked for generic algebraic fields. For 
many questions the real field R and all the p-adic fields Q, are on an equal footing 


conceptually, yet the Q, is often easier to handle. For example, one of the applications 
of the p-adic field Q, is to prove the nonexistence of solutions for certain Diophantine 
equations. 


5.2 Bruhat-Tits tree 


The Bruhat-Tits tree [59] is a special case of “building”, first introduced by Tits 
to study semi-simple Lie groups (in particular exceptional groups) using (incidence) 
geometry, and has since developed into a very powerful tool connecting group theory 
over arbitrary fields and geometry.'! (The original construction is in terms of simpli- 
cial complex (satisfying a set of axioms); however we follow the simpler construction 
in terms of lattices.) 


From the tensor network viewpoint, the motivation for studying the Bruhat-Tits 
tree is that it has the same tree structure as the one arises in the regular tessellation, 
but can furnish the full conformal group SL(2,Q,), much bigger than a discrete 
subgroup of SL(2, R), preserved by the tessellation geometrically embedded in the 


continuous AdS space. 

From the p-adic AdS/CFT, the Bruhat-Tits tree plays the role of the upper half 
plane whose boundary is the p-adic line Q, [60]. In this subsection we summarize 
the construction of the Bruhat-Tits tree, in particular motivating it from its role as 
the bulk of Q,, and prepare for the discussion on the SL(2, Qp) action on the tree in 
the next subsection. 


5.2.1 Bruhat-Tits tree as p-adic upper half plane 


The real field R is the boundary of the upper half plane H = SL(2,R)/SO(2,R). 
With coordinates H = {z = x +iy|x € R,y € R+}, it has SL(2, R)-invariant metric 


‘For good textbooks on buildings see e.g. [61-63]. 
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ds? = z (da? + dy”). An SL(2,R) action on a point z = x + iy on H would induce 
the same action on its boundary point x: 


az +b ax+b , & 
—> ==>. with 


L(2,R 12 
cz+d e std 4) Este, ) ine) 


If we replace the boundary space R by the p-adic field Q,, what would be its bulk, 
i.e. what is the p-adic version of the upper half plane? 


The analogy with the relation between H and its boundary R, shown in the 


following table 


upper half plane H | Bruhat-Tits tree H, 
Isometry group G SL(2, R) PGL(2, Q,) 
Isotopy group K SO(2, R) PGL(2, Z,) 
Boundary R Q 


suggests that one should replace R by Q, in the definition of H to give’? 


-~ .. PGL(2,Q,) 


We immediately see the difference from the real case. As the maximal compact 
subgroup of the isometry group PGL(2,Q,), PGL(2,Z,) is both open and closed 
(“clopen”) in PGL(2,Q,). Therefore, although Q, is a continuum, its bulk H, is 
actually discrete: it turns out that the p-adic upper half plane H, has the topology 
of an infinite (p+ 1)-valence tree (called Bruhat-Tits tree), instead of the continuous 


Q, ® Qp. This is in sharp contrast with the real case, in which both the bulk H and 
13 


its boundary R are continuous. 


5.2.2 Lattice construction of Bruhat-Tits tree 


Since H, has a discrete topology, we cannot simply give H, the coordinate zp = 
Lp + iYp in which zp E Qp and yp E Qp+ and write down the PGL(2, Qp) invariant 
metric on it. However, the coset expression (5.13) suggests that one can construct it 
in terms of the equivalence classes of lattices in Q, ® Qp. 


A lattice in Q, ® Q, generated by basis vectors ( f, g) 


i fi ~ [%1 . =o 
f= G g= a with fi gi € Qp (5.14) 


12The reason that we write PGL(2,Q,) instead of SL(2,Q,) is that the determinant of the 
SL(2, Qp) doesn’t have to be one, but only = 1 mod p: det g = p” + 1 with g € sl(2,Q,) 
13The Bruhat-Tits tree is an example of a Bethe lattice (the discrete lattice whose boundary is 


continuum), which has very interesting graph theory [64]. 
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which satisfies f 4 cg for Vc € Q, is defined as! 


(f,9) = {af +bgG|a,b€ Zp} (5.15) 


A PGL(2,Q,) transformation acts on the lattice via 


ER> Ph with q= (C n) E PGL(2,Q,) (5.16) 


The numerator PGL(2,Q,) in (5.13) acts transitively on the space of lattices in 
Qp ® Qp, and the denominator PGL(2, Zp) is the stabilizer of the integer lattice 
Zp Q Zp. Therefore H, defined as coset (5.13) is identical to the set of equivalence 
classes {A} of lattices in Q, ® Qp, where two lattices are equivalent, FD EP 
iff 


(f.7)=(-f,0-9) with TePGLOZ) (5.17) 


Note that since we use PGL instead of SL, we have 


(FD ~ MF) with EQ: (5.18) 


We denote the equivalence class of lattice generated by (f,g) as (f, g)). 


Now we are ready to define the metric on Hp, i.e. the distance between two points 


in H,, ie. two equivalence classes of lattice A and A’ in Q, ®Q,. First, A and A’ 
fis 


are defined as incident (i.e. directly connected) if and only i 
pACA'CA (5.19) 


We connect such a pair of nodes by an edge, which has distance d(A, A’) = 1. One 
can immediately check that each node has exactly (p+1) nearest neighbors, i.e. that 
H, has the topology of an infinite (p + 1)-valent tree (shown in Fig. 7 for p = 2), 
instead of a Q, ®Q, continuum. Then the distance between any two points is defined 


as the number of edges connecting them, which is invariant under PGL(2, Q,). 
For the real case, it is easy to visualize the relation between the upper half plane 
H and its boundary R. In the coordinate z = «+7 y (with which the hyperbolic metric 


is ds? = z (dz? + dy*)), taking lim,,9 gives the boundary point x and increasing y 
moves deeper in the bulk. An SL(2,R) action on a point z = z + iy on H would 


induce the same action on its boundary point x, see (5.12). In particular, in the study 
of holography, the constant y = € < 1 line is taken to be the cut-off surface. (The 


l4Here we see an important difference between the lattice on R @ R and the p-adic lattice. Since 
a,b € Zp (the unit ball inside Qp), which is continuous, the p-adic lattice forms a continuum, 
whereas the R @ R lattice is discrete. 

15To check that this definition is reflective, we note that since pA is equivalent to A, multiplying 
the equation above by p gives pA’ CAC A’. 
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choice of the cut-off surface is particularly crucial in the tensor network construction, 
see next Section.) What is the analogous situation for the p-adic case here? 


The lattice construction earlier only gives the topology of the tree. We would like 


to have a good coordinate system on the Bruhat-Tits tree H, such that its boundary 
is the p-adic field Q,, and the SL(2, Q,) action on a vertex in the bulk H, induces the 
same action on its boundary point. In particular, the natural cut-off surface from 


this coordinate system behaves well in the later tensor network construction. 


To fix a coordinate system in Hp, first choose the origin O. Since the tree is 


infinite and homogeneous, we can choose an arbitrary point as the origin O. We 
can then use the isomorphism PGL(2, Q,) to make O the equivalence class of integer 
lattice, i.e. its generators are 


O 


ae ; r 1 m 0 
(ogo) with f= (5) Jo = (o) (5.20) 
One first check that the (p + 1) neighbors of O are 


RGM CARON wih #20. eet G 


Note that these are precisely the neighbors that would arise when applying the Hecke 


T,(m) =m- l o) +m: J A (5.22) 


n=0 


operator 


on the node (( fo, go)). Now we use the projective equivalence (5.18) to fix the second 


1 
Applying Hecke operator iteratively then generates the entire tree, with all nodes 


1 
= 0 
vector g to go, i.e. rewrite the first node in (5.21) into (( i : ( J) 


having the form 


EEN Ear mer 62 


n=—N 


where we have used the projective equivalence (5.18) to fix g = 1. Figure 8 shows 
the Bruhat-Tits tree for p = 2, with the coordinate system (5.23). Note that since 
x truncates at p™, we can think of p™ as giving the accuracy level of a p-adic 
number 2"), i.e. the node (5.23) represents the equivalence class r™ + Lis: 


5.2.3 SL(2,Q,) action on Bruhat-Tits tree and the cut-off surface 


The coordinate system (5.23) assigns a vertex on the Bruhat-Tits tree two numbers 
p” and «™), In order to determine the radial direction and the boundary direction, 
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z E Qp 
oo 


(4f, g) p? 
(2f, 9) f g+ f p! 
(f, g) p? 

(F, 29) 


Figure 8. Bruhat-Tits tree for p = 2. Here f = P and g = go. 


and how to choose the cutoff surface in a manner suitable to holography (and anal- 
ogous to the upper half plane in the real case), we now study their behavior under 
the bulk PGL(2, Qp) action. 


A PGL(2, Q,) transformation acts on the lattice via 


FAO frd with y= (23) EPLO) — 6.24 


Given a vertex on the Bruhat-Tits tree with coordinate (5.23), under a PGL(2, Q,) 
action, it transforms as 


a(R), GON m () | (258) eon 


m| (ca +d)? 
ad — bc 


where 


p =p 


(5.26) 


P 


Recall that the starting point is « = ye an p” , with accuracy only up to 
level p”. The image is another bulk point at 


m'—1 


(m) l 
= a po bn p” with accuracy p” (5.27) 
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Now, as we go to the boundary, both m and m’ —> ov, 
m-1 lo) 
r™ = LD anp” — r= ` an p” (5.28) 
n=—N n=—N 


we have the expected boundary PGL(2, Q,) transformation: 


OO oO 


ax+b 
7 pz P cx+d Ps P ( ) 


Let’s compare to the real case. Under an SL(2,R) (5.12) on a point z = x +iy 


in the upper half plane H, 
x and y transforms as 


(ax+b)\(cx+d)+acy? yoo axt+b 


a aaa oe (5.30) 
yo 
y (cx +d} + (cy)? (cx +d)? 


Comparing this with the transformation of x”) in (5.29) and p™ in (5.26), we con- 
clude that p™ should play the role of the holographic direction, the analogue of y 
in the real case, whereas x‘ the role of the boundary direction. This also means 
that the cut-off surface should be a line of constant p”. We postpone the discussion 
on the implication of this choice for the tensor network construction until the next 
Section. 


5.2.4 Branches in Bruhat-Tits tree and the p-adic expansion 


Before studying in more detail the SL(2,Q,) action on the Bruhat-Tits tree, in 
particular the scaling transformation, we mention another way of visualizing the 
tree and understanding why it acts as the “holographic bulk” of the p-adic field Q,. 
First note that the p-adic expansion already has a natural tree structure. Consider 
a generic p-adic number 


boundary Q, : a= Ds anp with an € Fp (5.31) 
n=—N 
which can also be written like 
. 03 Ag Q1 Ag .A_1 A_gA_3 ...a_n (5.32) 


First, start from level p?, there are p choices for the coefficient ag, draw a node for 
each choice. The node with aj = 0 is then the origin O. Starting from each node 
(labeled by ag) at level p°, there are again p choices for a; — draw these p nodes at 
level p! and connect them to the node they start from. Moving up this way (and 
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also connecting all nodes at p? backwards to the node the correspond to Op~! + 0), 
one draws an infinite (p + 1)-valent tree starting from level p™t. Moving backwards 
towards negative power p`” then give the entire Bruhat-Tits tree. At level m, a node 
correspond to 


m-1 
Bulk tree at level m : r™ = ` an p” with an € F, (5.33) 
n=—N 


Then for a given p-adic number, starting from the lowest power and at each level 
following the branch corresponding to a, in the expansion (5.31), we have a one- 
to-one map between a p-adic number and a branch in this (p + 1)-valent tree. To 
summarize, each node in the bulk Bruhat-Tits tree has two label: (z = x'™, z = p™), 
where p™ gives the accuracy level and x”) a p-adic number to the accuracy p”™, i.e. 
it represents the equivalence class x”) + pL: 

T reEQ@ 


pth 0  _ 
| 


Figure 9. Bruhat-Tits tree for p = 2. The dilatation D on the tree slides the branches 
along the main brunch. 


The non-zero elements in Q, can be grouped according to the leading term in 
the p-adic expansion: 


Q; = Jre U, : tp=p u with vez and uecU, (5.34) 


neZ 


— 36 — 


ie. u = > oan p” with ao # 0. The set p"U, for each n € Z form a subtree with 
the root at: 


Points on the main branch: a” =p”.-0 with accuracy p” (5.35) 


The line connecting all the re is then the main branch, running from n — —oo to 
n — +00. 


5.2.5 Scaling transformation and scaling dimension 


It is quite remarkable that the Bruhat-Tits tree, though discrete, can furnish the 
full conformal group PGL(2,Q,). This allows us to study the function on the tree 
which has a definite quantum number. Given the Iwasawa decomposition G = NAK, 
where N is the Borel subgroup, A the dilation, and K maximal compact subgroup 
SL(2, Zp), it is enough to look at their actions separately. Since we will not use the 
detailed action later, we leave this to the Appendix, and only discuss the dilatation 
by p” here. 


Under a dilatation by p”, a vertex on the tree transforms as 


p= (Pr eal ECT «PE. Pa py 38 


The action moves the branches along the main branch, shown in Figure 9. 


5.3 p-adic Fourier transform and p-adic wavelet transform 


In this subsection we review some basis on p-adic integration, p-adic Fourier trans- 
form, and p-adic wavelet transform, which will be needed later when we discuss bulk 
operator reconstruction in p-adic AdS/CF'T in section 6.3. For more on p-adic anal- 
ysis see the textbook [48]. The summary on p-adic integration and Fourier transform 
here follows mainly the chapter 3 of [65]. 

5.3.1 p-adic integration 


The integration measure is not unique. There are several commonly used measures 
for the p-adic integration. First, the additive measure dx over the Q, line is defined 
by demanding the following translation and scaling behavior 


d(x +a) = dz and d(ax) = |a|, dx (5.37) 


Usually it is normalized over the “unit ball” (defined in (5.11)) 


[ dz =1 (5.38) 
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Not surprisingly, the integration over the entire Qp using the additive measure dx 
diverges. The multiplicative measure d*x over Q% can be defined by 


p dz 
d* (x) = —~— 5.39 
(=) p—1|z|p l ) 
which satisfies d* (ax) = d*a and fņxd%x = 1. Finally, the Patterson-Sullivan 
measure d(x) over P+(Qp) is defined as 


eae ma (5.40) 
ae otherwise 


The distance as seen by this measure agrees with the one computed by counting 
steps in the BT tree. 


5.3.2 p-adic Fourier transform 


To define the p-adic Fourier transformation, we need the analogue of e?7"**, i.e. the 
additive characters on Q,. This is defined as'® 


wose "m (5.41) 


where [|x] is the “fractional part” of the p-adic number z, defined as 


[x] = 5 anp” (5.42) 


and |x] = 0 for x € Z,, i.e. |x|, < 1 (hence the name “fractional part”). It is easy 
to check that x(£)Xr(y) = xk(£ + y) and |x;(x)|, = 1. Integrating the additive 
character X(x) over the unit ball Z, gives the characteristic function y(k) of Zp in 


Qp: 
, 1 keZ 
I dx e7? ka = y(k) = i p (5.43) 
Zp 


0 otherwise 


Note that both the additive character x, and the characteristic function y(k) of Zp 
in Q, are locally constant functions.” 

For the study of p-adic Fourier transform, we summarize a few useful integrals, 
which are interesting in their own right. First, the integration over bigger balls p” Zp 


gives the characteristic function of p"Z, in Qp, i.e. y(p"k): 


1 
L dre lkal = PH (5.44) 
p” Lp 


16 The sign difference in the exponents between the real character and the p-adic ones is important 
for constructing the adelic product. 
17 As a contrast, recall that for a function from R to R, “locally constant” implies “constant”. 
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From this we can compute the integration over “shells” 


on = |Klp < po? 
/ a dz e7? [ke] — — |Klp = pet (5.45) 
O? 0 ker 


where U, is the “unit sphere” of p-adic numbers defined in (5.10). This in turn gives 
i dx e’ ke] — _-() (5.46) 
Qp\Zp 


The p-adic functions that allow a well-defined (i.e. invertible) Fourier transforms 
are locally constant (i.e. constant within each p"Z*) functions with compact support 
or sufficiently fast asymptotic decay. The p-adic Fourier transform and its inverse is 
then 

fk) = I wre T and f(x) = f dk f(k)enl (5.47) 
P P 

The Fourier transform of the following two functions will be needed later. The 

characteristic function y(x) of Zp in Q, is invariant under Fourier transform, i.e. 


lz) = qlz) (5.48) 


hence y(x) is also called p-adic Gaussian. 


5.3.3 p-adic wavelet transform 

Let’s focus on the one-dimensional case. Given a p-adic mother wavelet y(x), which 
for our purpose is taken to be a complex function with p-adic argument x, her p-adic 
daughter wavelets can be defined in analogue to the continuous case (3.13) 


Wa,s(Z) = 


1 r—-a 
Y 
y lslp l 2 


This set of daughter wavelets forms the basis for the p-adic wavelet transform. Then 


) z,a,s E Qp (5.49) 


the p-adic wavelet transform of a function f(x) with z € Qp using this set of the 


wavelet basis is i 
r—a 
dx wr )f(z) (5.50) 
Qp V [slp . 


whose inverse transform is given by 


W;(a,s) = 


1 ds 
j2)= C, L; isp i. da W (a, s) a,s(X) (5.51) 
with TE 
o= f a OF sa 


P 
dsda 
[s12 
p-adic affine group (whose action is x —> sx +a). For the inverse transform to exist, 


Note that the integration measure is precisely the left-invariant measure of the 


we need the same “admissibility condition” Cy < +00, as in the real case. 
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6 Tree-type tensor network and p-adic AdS/CFT 


In the previous section, we have described how tensor networks naturally re-enact 
various features of the AdS/CFT dictionary. These include the emergence of the 
HKLL relation in reconstructing bulk operators from boundary operators, and the 
relationship between global symmetries in the boundary and gauge symmetries in 
the bulk. Having discussed also the symmetries of the Bruhat-Tits tree, here we 
would like to see how the p-adic AdS/CFT developed in [26, 29] re-emerges in the 
tensor network. 


6.1 Emergence of p-adic AdS/CFT 


In [26, 29] the authors generalized the AdS/CFT correspondence to a BT-tree/p- 
adic theory correspondence. In this case, it is not clear what the boundary theory 
is. For simplicity, suppose we focus on BT trees defined on a two dimensional graph. 
Then given the bulk isometry of the BT tree, one can be assured via the usual 
understanding of the AdS/CFT dictionary that the boundary theory is one dimen- 
sional. Moreover, the theory resides in the space Q,,, and it should be invariant under 


SL(2, Qp). 


Since the correspondence is guessed based on the AdS/CFT correspondence, not 
every feature has an immediate counterpart. Indeed, one important ingredient in the 
AdS theory, namely gravitational excitations, does not have an obvious counterpart 
in the BT tree. There is an attempt to introduce a discrete version of Einstein theory 
in the BT tree [30]. However, in the following, we will mainly focus on the scalars, 
which have a more straightforward generalization on the discrete graph. 

In [26, 29], it is assumed that the action describing a bulk scalar field assumes 
the form 
p(w) d*(v) 


3l + 14 4! pre , (6.1) 


I= Y -2o + m3)o(v) + m 


v 


where v denotes vertices on the tree, and O is the graph Laplacian defined in (3.7). 


The setup is thus very similar to the continuous case, except that now one proceeds 
with solving for modes on the graph, and obtains propagators and mode expansion 
of the fields. 

As we have seen in the discussion in section 3, independently of the precise 
graph the tensor network lives in, relationships between bulk and boundary operators 
acquire a relationship that is based naturally on solutions of the graph Laplacian. 
Here, therefore, as we reconstruct bulk operators from boundary operators according 
to the HKLL formula now transplanted to the BT tree, the expressions that we 
obtain, namely (3.5) becomes exactly what one would have recovered according to 
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the BT tree/p-adic CFT correspondence. The relationship between mass and scaling 
dimension is also exactly as given in [26, 29]. 


There is one potential confusion concerning the signature of the theory. In the 
usual discussion of the tensor network, such as the MERA and the MPSs, they are 
descriptions of wavefunctions, which are defined on a given time slice. In the p- 
adic version of the AdS/CFT [26, 29], the construction is a generalization of the 
Euclidean version of the AdS/CFT. How then should the tensor network be un- 
derstood in this case? To understand what happens, we recall works in the tensor 
network literature where classical partition functions of statistical models admit a 
tensor network description, see for example [66, 67]. These tensor networks repre- 
sentations of partition functions, which are equivalent to Euclidean path-integrals of 
the quantum model in one higher dimension, can also be coarse grained, which are 
linear maps of the constituent tensors. If every step of the coarse graining is kept 
so that those linear maps form layers of tensors in an extra dimension, we obtain a 
tensor network in one extra dimension. In other words, there is no mystery about 
the Euclidean version of the tensor network/geometry correspondence. It would have 
the same form as MERA except that what was previously physical dangling legs are 
now contracted among themselves in the Euclidean path-integral. The bulk operator 
boundary operator reconstruction is clearly independent of whether the boundary 
layers of legs are contracted among themselves, and therefore do not get in the way 
of our understanding. 

The same analysis regarding wavelet transform that applies in the usual AdS/CFT 
correspondence would also resurface here, now as a direct result of the tensor net- 
work construction rather than part of the conjectured dictionary in the AdS/CFT 
correspondence. This is discussed in the following. 


6.2 Scaling symmetry — cutoff surface and dynamics 


There is something missing in the discussion of global symmetries and tensor net- 
works in subsection 4.1. A conformal field theory has a family of primary operators 
whose conformal dimensions are determined by the dynamics of the theory. The 
states are arranged into representations according to the primaries present in the 
theory. Therefore it is not possible that the tensor network realizes the conformal 
symmetry purely as a geometric symmetry, or we would come to the same conclusion 
as we did with translation symmetry, with no input that derives from the dynamics. 


On the other hand, it is already understood where the information of conformal 
dimensions are hidden. They are precisely encoded in the values of the tensor ele- 
ments. In the case of MERA type tensor networks, where moving through layers of 
tensors realizes coarse graining, the procedure of recovering conformal dimensions of 
primary operators in MERA is almost exactly identical to the current holographic 


= Al 


code. At a fixed point, the operator acting on the legs of a given layer is mapped by 
the tensor to a “coarse-grained” operator. At a fixed point for a primary operator 
however, it is rescaled under coarse graining, giving 


T [is] [js] =, s+1 
2, Ta is(a) js(b) Ties = AO; is+1(a)fs+1? (6.2) 
lis],ljs] 


where \ = A~4, see Figure 10. The symbol [i,] denotes a set of indices of which only 
i,(a) is contracted with ©. Here A depends on the ratio of coarse graining where A 
sites are merged to one. 
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Figure 10. This is an illustration of equation (6.2). 


The operator maps that we discussed in (3.1) is closely related to the expression 
above. How should we fit together then that conformal symmetry is a geometric 
symmetry of the bulk network and yet contains dynamical information stored in the 
tensors? It turns out the key point lies in the choice of the cut-off surface, more 
precisely, in the choice of the holographic direction along which one goes to the 
boundary. In the usual AdS/CFT correspondence, the boundary limit often leads to 
divergence in the action, and it is a common procedure to regulate by introducing 
a cut-off surface. In many practical applications of the tensor network, the network 
involved is most commonly finite and ends somewhere. The issue is thus about what 
is a natural choice of the cut-off surface. 


In AdS» for example, SL(2, R) transformations act in the Poincare coordinates 
according to (5.30). To study the CFT on the real line such that the coordinate x 
transforms according to (5.12), it is thus natural to choose a constant y = e surface 


as the cut-off surface. As discussed in the previous section, the accuracy p™ plays 
the role of y, and x”) the analogue of x in the BT tree. Therefore it is natural to 
choose a constant accuracy p™ as a cutoff surface in the BT tree. 

One important feature on the cutoff surface is that the number of physical dan- 
gling links are not evenly distributed. Instead they get denser as one moves along 
the surface. The distance of the dangling legs from the main-branch increases as we 
move along the cutoff surface. Therefore, as we apply a rescaling transformation at 
the boundary, we would hop from branches closer to the main branch to branches 
further away from the main branch. Under such a rescaling therefore, operators 
would have to be sandwiched in the same manner as in (6.2) and thus recover its 
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scaling dimension. As illustrated in Figure 11, an operator that acts on link 2 would 
be transformed to an operator that acts on link 4, which is related to the original 
operator as if it is sandwiched between an extra layer of tensor (see the box inside 
Figure 11. Therefore in the tensor network reconstruction, the cutoff surface is key 
to recovering the dynamical information that comes with scaling transformation. 


Figure 11. We illustrate the analogue of a Poincaré cutoff surface in the BT tree. Tensors 
move down the tree under scaling transformation. We indicate with green arrows a few 
examples how the links move. Such a transformation corresponds to transformation of 
operators such that it is sandwiched between a network tensor. This is depicted graphically 
in the box. 


6.3 p-adic wavelet transformation and bulk operator reconstruction 


Having fixed a natural orientation in the tensor network, in this sub-section we extend 
the arguments of section (3.1.4) to the p-adic case, and derive the p-adic version of 
the bulk operator reconstruction and its connection to p-adic wavelet transform. 


Recall that in the real case, the “smearing function” (3.26) used for the bulk op- 
erator reconstruction is related to the bulk-to-boundary propagator by first replacing 
the conformal dimension A in the exponent by d — A and then multiplying the © 
function in order to restrict the integration to the causal patch. In this subsection, 
we will see that the p-adic bulk operator reconstruction follows exactly the same 
pattern. Moreover, just as in the real case, this bulk operator can be regarded as 
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the p-adic wavelet transform of the boundary operator. (In this subsection we will 
mostly focus on the case of d = 1.) 


In the p-adic case, the bulk-to-boundary propagator of a scalar with mass m is 
[29] 


o lzlp 7 
Deal (oo Roar) a 


We remind that here x and y are coordinates of the boundary Q, line, and z = p” 
with m € Z denotes the holographic direction in the BT tree, and gives the accuracy 
of the p-adic expansion. In particular, |z|, = p™™, and |z|, —> 0 as one approaches 
the boundary. The supremum norm sup{|z|,,|% — yla} = |2|p if |z|p > |z — y|» and 
|x — y|, otherwise. In equation (6.3) the conformal dimension A is related to the 


. 2 1 . . 
mass via m* = — and M, is a numeric constant. 
&@(A-1)& (4)? p 


In addition to replacing the exponent A in the p-adic bulk-to-boundary propa- 
gator (6.3), we multiply the propagator by the analogue of the © function in order 
to restricts the integration in the manner as in (3.26). This gives the p-adic analogue 
of the “smearing function” 


Kala, zly) = Np lela AE) (6.4) 


where y is the characteristic function of Z, in Q, defined in equation (5.43). 


The linear term of the bulk operator reconstruction, using the “smearing func- 
tion” (6.4), is then 
dlez) = | dy Kyl2, zly) OW) (6.5) 
P 
Just as in the real case, we can interpret (6.5) as a p-adic wavelet transform (5.50). 
However, the main difference from the real case is that now the choice of the mother 
wavelet does not depend on the conformal dimension A of the boundary operator. 
For any operator, there is a “universal” mother wavelet 


(a) = Np (2) (6.6) 


The conformal dimension A only enters through the relation between her daughter 
wavelets and the “smearing function”: 


K,(2,2|y) = lel vt. (y) (6.7) 


which in turn gives the mapping between the wavelet transform Wo(z,z) and the 
bulk field: 
At 
pplz, 2) = |z| “Wola, z) (6.8) 
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Similar to the real case in (3.31), our mother wavelet (6.6) does not obey the 
standard admissibility condition since 


2 dk 2 2 f E 
Co NG EWP ANG fg te (6.9) 


where we have used the fact that y is self-dual under Fourier transform. Again, this 
does not pose a problem for us, because the physical inverse transform we use to 
obtain the boundary field O(y) is 


Jate | & oe ax (z, zly) (6.10) 


which has an extra dressing factor of |z|?4~!, due to (6.7) and (6.8). We will now 
show that this can make the normalization C finite and the inverse transform (6.10) 
valid. Note that similar to the real case (3.26), the function y restricts the integration 
within the causal patch. (In the Bruhat-Tits tree, the causal future of a node (z, z) 
is defined as the subbranch of the tree rooted at this node.!® ) 

First we compute the Fourier transform of the p-adic “smearing” function (6.4), 
using the additive character (5.41): 


Kalm aly) =N lz | dhe 6-0 (kz) (6.11) 
p 
Note that since the function y is self-dual under the Fourier transform, in the mo- 
mentum domain the restriction to time-like momenta is also implemented by the y 
function. Plugging (6.11) into the inverse transform (6.10) we get 


NRG) | ax'Gyla! = YOO) (6.12) 
with 
G,(2' — y) = | dk e eE (6.13) 


P 


the p-adic analogue of the kernel in (3.37), and!° 


P E = p=) 
G= | dze? te)? = f deep = forv > 0 6.14) 
p 3 | |p ( ) Zp | lo pl — p-?”) ( 
with v = A — z. Therefore, for all v > 0, i.e. for massive scalars, the p-adic 


normalization C, is finite and the inverse transform is valid. 


18The causal structure on the Bruhat-Tits tree was first studied in [68], in the context of hierar- 
chical structure of the inflation. 
19Note that the integral in Cp only converges when v > 0. 
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We have thus shown that the p-adic bulk reconstruction works in the same 
manner as the real case: the bulk reconstruction can be realized as a continuous 
wavelet transform, whereas the inverse, i.e. obtaining the boundary field from the 
bulk one, is related to the standard inverse wavelet transform by an extra dressing 
factor of 22” 


There is one important difference between the real and the p-adic cases. In the 
derivation of the bulk reconstruction via a p-adic wavelet transform, we have used 
the continuous wavelet transform. However, in essence, the p-adic wavelet transform 
is discrete.” The reason is the following. First of all, even though we started with a 
continuous wavelet transform in (5.50), the resulting scaling parameter is |z|,, instead 
of z. And since the integration measure in the inverse transform (5.51) is invariant 
under the affine group, the integration is actually equivalent to the sum over the 


Bruhat-Tits tree: 
fe A dx f(x, z) = 5D, Fe) (6.15) 
|z| aa 


er 
Recall that the Bruhat-Tits tree can be viewed as the discrete, holographic, bulk 


whose boundary is the continuous line of Q,. The emergence of the discrete p- 
adic wavelet transform starting from a continuous one can be viewed as the mirror 
statement of the above. This thus lends strong support to the connection between 
the bulk reconstruction and the wavelet transform, with the identification between 
the coordinate in the holographic direction in the bulk reconstruction and the scaling 
parameter in the wavelet transform. 


Finally, we also note that the bulk reconstruction discussed in [26] is the Eu- 
clidean version of recovering the bulk field for specified Dirichlet boundary condition, 
which use the non-normalizable bulk-boundary propagator, hence it is a different sit- 
uation from the HKLL relation. 


6.4 Correlation functions 


We have postponed the discussion of correlation functions which was described in 
detail in [26, 29]. Correlation functions between local operators in the boundary is ob- 
tained. They are built up from tree graphs involving the bulk-boundary propagators 
and bulk-bulk propagators. These propagators, despite the need for an appropriate 
choice of normalization, are indeed proportional simply to pô»), where d(u, v) is 
the graph distance between the vertices u,v, which is unique on a tree, and A is the 
analogue of a conformal weight in the SL(2,Q,) invariant theory at the boundary. 
Let us pause to make the definition more precise. 


20The coincidence of the continuous and discrete p-adic wavelet transform was already shown for 
a particular type of wavelets in [78]. 
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A primary field O4 (x) of PSL(2,Q,) with conformal weight A is defined as [73] 


ax+b 
O' = 
oe + a ( 
Now take for example the two point function. The functional form is completely 


fixed by the PSL(2,Q,) symmetry. Following the standard AdS/CFT dictionary, it 
gives [26, 29] 


det y 
cr+d 3 


A 
Oa (z) where y= C € PGL(2,Q,) (6.16) 


C 


|z1 — a 


(Oa(21)Oa(22)) = (6.17) 


where ||, is the p-adic norm. Up to normalization that removes the divergence as the 
cutoff surface is taken off to infinite accuracy m — oo, this is basically evaluating 
pruer22) 


£1, £2 through the BT tree. This can be seen using Figure 9 of [26]. This result is 


where d(x1, £2) is the graph distance between the two boundary points 


derived via the AdS/CFT dictionary, and we would like to see how this could emerge 
from the tensor network without a bulk action or the AdS/CFT dictionary that 
we could invoke. Computing the correlation function is somewhat more complicated 
than recovering the HKLL relation. The reason is that the computation of correlation 
functions involves insertion of external operators at the boundary physical legs. A 
priori, it is not clear how to push the action of a generic operator into the interior. 

There are some special cases in which this simplifies so that an explicit calcula- 
tion is possible. For example, correlation functions can be computed approximately 
when the tensors are close to perfection [21]. Since it appears that the dimension 
of an operator can be recovered in the current discussion in a manner completely 
analogous to MERA according to the explanation near Figure 11 and (6.2), it is 
possible to compute the correlation functions in an analogous manner as in MERA. 
One consistency condition following from physical expectation is that the identity 
operator should scale trivially under (6.2). If such is the case, we have (see Figure 
12) 


lis] ts 
dT a = OE EE (6.18) 


Figure 12. This illustrates the consistency condition in equation (6.18). 
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Now using both (6.2) and (6.18), one would recover a two-point correlation func- 
tion given by 
(Oa(i1)Oa(iz)) = AA), (6.19) 


where 
A= D Taia Tj sim Oi aja (2) Orn Bhim (6) (6.20) 
[im], [3m] ,Jm—1 
see Figure 13. Equation (6.19) would thus recover the correct behaviour (6.17). The 
computation of the two-point function is illustrated in Figure 14. These considera- 
tions can be readily generalized to computing three point functions, see [77]. 


O 


Figure 13. This is an illustration of the contraction in equation (6.19). 


Having got this far, one is tempted to believe that an interacting scalar theory 
living in the discrete graph has emerged. But this theory is not yet one that is 
describable by (6.1). There are some extra consistency conditions clearly required to 
make this possible. 

First of all, if the propagator following from (3.5) is to have the same scalar 
mass consistent with the equation (3.10) independently of the direction the particle is 
propagating, then it is necessary that the legs of each node are permutation invariant. 
One implication is that combining with (6.18) and taking p = 2, the tensors satisfy 

T* 


abc 


Tabe = Dech; (6.21) 


and this is true for any choice of two of the three legs to be contracted in the above 
equation. We note that this is the closest analogue to the perfect tensor condition 
for a tensor with an odd number of legs. One might worry that the p = 2 BT 
network would fail to recover any connected correlation function, for the same reasons 
discussed in [21]. Miraculously however, precisely for p = 2 on the tree one can see 
that the argument of [21] breaks down, and the connected correlation function is still 
generically non-vanishing without a need to perturb them away from perfection. For 
p > 2 the consistency condition gets less stringent in comparison. 
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On top of that, we note that in equation (6.2), some conformal dimension is 
obtained. However, information about primary operators have already showed up 
elsewhere: namely in the HKLL relation. There are two consistency conditions that 
follow from the usual AdS/CFT dictionary. First the basis of primary operators 
that is obtained in the scaling relation (6.2) should coincide with the basis that 
diagonalizes the matrix a‘ defined in (3.1). Secondly, what we have denoted ø in 
(3.4) should be related to A by 

A+o=d. (6.22) 


At first sight this is not necessarily true. However, if we adopt a “symmetric re- 
construction”, in which operators acting on a given leg of a p+ 1 valent is pushed 
to the other p-legs, we can prove that in the weakly coupled limit (where Og. Bp 
controlling interaction is much smaller than the propagator term, the basis vectors 
for operator pushing and that of the scaling relation coincides, and that (6.22) is 
satisfied automatically. These results will be reported in a publication currently in 
preparation. 
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Figure 14. This illustrates the two-point correlation function on the 2-adic tree. 


7 Summary and Discussion 


In this paper, we examined the tensor network/geometry correspondence in further 
detail. Intent on finding a counterpart of the HKLL relation, we find that such rela- 


— 49 — 


tions do exist in the tensor network, mirroring closely the AdS/CFT dictionary. In 
particular, the smearing function in a homogenous network readily admits an inter- 
pretation as a solution of the graph Klein-Gordon equation, with masses determined 
by the specific matrix elements of the tensors. Moreover, such relations naturally 
give rise to interaction vertices that supply the building blocks for loop corrections, 
mirroring the systematics of 1/N corrections in the AdS bulk. 

The story becomes richer when there are more symmetries under control. To 
that end, we turn our attention to the recent p-adic AdS/CFT correspondence. By 
putting our tensors in the BT tree, we find that the general properties of tensor 
networks pass through naturally to the BT tree, but now with the advantage that 
there is a full SL(2,Q,) symmetry. With that, we could explicitly construct these 
smearing functions in terms of Green’s functions of graph Laplacians, and identify 
the scaling transformation that allows the scaling dimension of primary operators 
of the p-adic CFT described by the tensor network to be read off. The systematics 
generically re-enact the proposed dictionary of the p-adic AdS/CFT correspondence. 
Let us emphasize here as we did in the main text, that we could interpret the tensor 
network as a wavefunction at a given time slice, or as a path-integral of a Euclidean 
(Lorentzian) partition function by contracting the external legs with corresponding 
matrices that are the tensor network approximation of e7™” (e*Ħ). The conclusions 
reached about generic tensor networks are insensitive to the treatment at the bound- 
ary. 

The current work is only a baby step in elucidating the relation between tensor 
networks and the AdS/CFT correspondence, if the relation is indeed more than a 
pictorial analogy. We list some of the imminent questions and generalizations in the 
following. 


7.1 Implication to p-adic AdS/CFT 


One reason that p-adic AdS/CFT is much less developed than the real version is that 
there hasn’t been a well-defined, but non-trivial, “p-adic CFT”. One of the main 
difficulties is that the ordinary derivative ð, of the real case does not exist for the 
p-adic field Q,. To proceed, there are two alternatives. One is to use the “Vladimirov 


derivative” : f(y) — f(x) 
Dyf (2) = f ay 


which defines the derivative via integration, and hence can be regarded as the p-adic 


(7.1) 


analogue of Cauchy’s contour integral in complex analysis. The resulting field theory 
with a Lagrangian based on the “Vladimirov derivative” is therefore highly non-local. 
Until recently, only the free massless scalar has been studied [69-71]; see, however, 
the recent work [72] for an interacting example. 

The other choice is to simply forgo the use of derivative, as in [73]. The “p-adic 
CFT” defined this way would have only primary fields (defined using global SL(2, Q,) 
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as in (6.16)), and no descendants. The unitary condition is then automatically 
satisfied and therefore cannot impose any constraints on the conformal weights of 
allowed primaries. In fact, the form of all correlation functions are fixed and the 
crossing symmetry is automatically satisfied, therefore it also cannot impose any 
constraints on the allowed primaries.” 


With the result of the present paper, we can regard the tensor network living 
on the Bruhat-Tits tree as a concrete realization of a “p-adic CFT”, with different 
choices of the tensors corresponding to different matter content. Since a priori, any 
conformal dimension is allowed, and we do not need the derivative or Lagrangian to 
define the theory, this is more in line with the second, more algebraic, approach. 


This then allows us to go beyond the free massless scalar. The next simplest 
examples would be Wy, minimal models,”” which also has the benefit that it is 
dual, by a weak/weak AdS/CFT correspondence, to a bulk higher spin gravity [74]. 
The existence of a weakly-coupled bulk dual living on the Bruhat-Tits tree then 
imposes additional constraints on the boundary “p-adic CFT”. This is currently 
under investigation. 


On the other hand, it was observed in the continuous version of the AdS/CFT 
correspondence that the bulk dual of the (global) conformal blocks are geodesic Wit- 
ten diagrams [75, 76]. It would be interesting to look for the analogous statement 
in p-adic AdS/CFT, which would allow us to read off the p-adic CFT data directly 
from the values of the residing tensors. We hope to report on this soon. 


7.2 Further explorations in tensor networks 


In this paper, we have studied very specific types of graphs. Namely they are un- 
weighted graphs with a fixed valency. Our language should admit generalizations 
to these more general types of graphs. Let us comment briefly on two interesting 
extensions. 


e For tree graphs with multiple valency, it can be described by a substitution 
matrix. A substitution matrix M is a N x N dimensional matrix, where N is 
the number of different species, and M;; is the number of descendants belonging 
to vertices of type 7 branching out from a vertex of type i . For example, if we 
have only two kind of vertices, we have 


M= & (79) 


21The modular invariance on higher genus Riemann surfaces might impose more constraints but 
this has not been studied. 
22Tensor network has been used to study the Ising model and three-state Potts model in the real 


case in [66]. 
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where the matrix elements are positive integers. One solution to the (un- 
weighted) graph Laplacian is 


Gola) = NA™K", (7.3) 


where M is some suitable normalization and m + n = d(a, O), which is the 
distance between a node a and the root of the graph O.” To satisfy the graph 
Laplacian equation with the same mass, we need to put in an extra constraint 
so that the function satisfies the same equation independently of the type of 
node the Laplacian acts on. The constraint is given by 


1 1 
ch-+ dk +5 —(e+d+1) =ah+bk+ = —(a+b+1). (7.4) 


It is not yet clear what such constraints could mean as far as the bulk interacting 
— theory is concerned. It should be studied systematically. 


© e Another generalization is to include weighted graphs. In this paper, we have 
i assumed that the network is a homogenous network and preserves as much 
symmetry as the graph. In general when we start discussing fluctuations around 
a background, then it seems natural that different edges should in general have 
different weights. A first step has been taken in [30], where the weights are 
interpreted as metric fluctuations. The Laplacian also depends on these weights 
by 


glv) = X Jwlolu) — (0), (7.5) 


UNV 
where Jw is the weight on the edge connecting the vertices u and v. It is an 
important question to understand how these degrees of freedom emerge from 


the tensor network. 


Apart from these generalizations, there are other important questions. In par- 
ticular, it is evident that to recover a bulk theory that can be interpreted as a local 
interacting theory, there are many constraints that need to be imposed on the tensors. 
It is necessary to study these constraints systematically, and ask for the minimal set 
of data needed to define the bulk theory. It would also be important to study the 
distribution of eigenvalues that control the conformal dimensions of operators in the 
large bond dimension limit.” Such a discussion is pertinent to understanding when 
a large gap in conformal dimension is emerging, which is believed to be an essential 
condition for the emergence of a semi-classical and local holographic dual [79]. 


Second, we need a theory of gravity in the bulk, and to that end, we need to 
find a systematic way to describe fields carrying spin. As far as gravity is concerned, 


?3Roughly speaking, a rooted tree contains a special point from which the tree grows outwards. 
24We thank Wei Song for the suggestion. 
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one way to identify fluctuations of these tensors with gravitational excitations is 
based on diffeomorphism. The idea is that the bulk graviton is supposedly dual 
to the boundary stress tensor. One could therefore ask the following question: for 
given local diffeomorphism at the boundary, such as a local translation that can 
be effected by inserting the (exponentiated) stress tensor ei J dre(x)Tox we can look 
for the corresponding transformation in the bulk that would have (approximately) 
the same effect. This would recover 6g, according to the JT), inserted. The bulk 
transformation in general corresponds to the exchange of edges (which is particularly 
apparent in the case of a tree). It thus fits with a metric carrying two indices. This 
is currently under investigation. 

Last but not least, as we already mentioned in section 2.2, recovering bulk dy- 
namics from boundary dynamics is an extremely important question. We hope to 
report on these issues again in a future publication. 
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Appendix 


A HKLL-like relations and fusion in a tensor network near 
perfection 


In this appendix, we will demonstrate how the discussion of HKLL-like relations and 
fusion matrices in section 3.2 are related to the discussion in [21]. To that end, like 
[21], we will have to work in the limit where every tensor is close to being a perfect 
tensor. Recall that a perfect tensor is one satisfying 


Tal = big, (A.1) 


where J and a each denote an arbitrary set of n indices among the 2n indices of T. 
Now we will consider tensors at each site taking the form 


T=T+et, (A.2) 


=5)3= 


where T is a perfect tensor, and t is some arbitrary tensor and e an infinitesimal 
parameter. 

First we would like to show that the coefficients a for the linear term can at 
most be of order O(e). 

We take a basis of operators acting on each link to be traceless apart from the 
identity matrix. (i.e. for bond dimension D = 2 we take Pauli matrices oz, oy, 0, as 
the basis. For general cases we take a generalized set of Pauli matrices generated by 
X,Z satisfying 

X.Z = exp(2ni/D)Z.X. (A.3) 

Now we write equation (3.1) explicitly as matrix multiplication. Wlog, we act 

OK on the first leg 


Oaai (T + El) piace = >, aF (1, DT > A E E + 
bilas 


af (1, l, lo)(T + Ch) onBi Biz OF OZ see (A.4) 


Big Qla 


For generic t, even though T ceases to be unitary, its inverse exists independently of 
how we split the legs into two groups. Therefore, in the equation (3.1) we are free to 
choose L/2 legs out of the L legs on which non-trivial operators could act on. 

Consider to zeroth order in e. We would like to extract the index a/‘(1,1) for 
some leg l. On the lhs, we contract all the legs except 1 and / with Tt. Clearly on 
the left we have OF ða- On the rhs, we have 


&œıyı 
Maf (1, Dary Ce + 
L/2-1 


K sk l fj 
OF, Tyg Ay C1) €25 67 jot) Lg Qel esseri da ye Bye joa I] Old + 
i 


L/2 


K t II L 
ar eiaa l E1, E2; +: €L/2) Toy ay er ezer jan Dyr y mepeg ezja Ocie (a0 
i 


where we have discarded any term containing trace trO. Note also that the expression 
above does not involve operators acting on more than half of the legs of T. That’s 
because T is perfect and we can always pick leg 1 together with another L/2— 1 legs 
as the input and the other L/2 legs as the output of the unitary map defined by T 
with this particular grouping of indices. Now if we trace over a,y, as well, the last 
two terms above vanish, along side the lhs in which we take trace of OX. Therefore 
we have to conclude that all a vanish to zeroth order in e. 

Now we can repeat the proof for a‘,(1, l, l2) by not tracing out two of the legs 
l, l2, but using the fact that a* itself vanishes. Again we would conclude that both 
the Ihs and all the terms in the rhs vanish except for a,(1,l,,l2), showing that it 
jy ae Therefore 
the first coefficient that would show up at e° has to be aj .. toe This is indeed the 


has to vanish. This argument can be iterated all the way to ap. 


case in every explicit example of the perfect tensor. 
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A.1 Connection with transfer matrices 


Here we would like to show that the transfer matrices we have defined in [21] are 
closely related to the a¥ that we have observed here. Return to equation (A.4). Now 
contract all the indices except link l with (7’+ et)’, and extract the order € term. We 


have 
K K I K 
Cee ee wie T) Oni = ay (1, DO; U + O Tah €l, €2,77? eLj2) 
L/2 
T L 
Ce eer er ree +te T) I] One: (A.6) 
i 


The lhs is precisely the transfer matrix contracted at one end to the operator 
OK. The rhs contains both the linear term and a O"/? term. The basis for conformal 
operator in this limit is a mixture of the basis defined by the linear term a‘ and its 
complement term containing O//, since to zeroth order in € they are interchangeable. 


A.2 Connection with fusion matrices 


Here we inspect also the meaning of the fusion matrix in this light. Again we start 
from (A.4) and contract with (T + et)" all the indices except two legs l,l. This 
gives, to order € the fusion matrix with one end contracted to O* on the right, and 
on the left we have 

ak (1, 11) O2,, br, Ooty Big + ly = ly + 


o (l,l l2)OL, p, O2,,0,, + 


Ql biz Aly 


K Ny Tig 
OT, Re ee ee li, l2, €1, €2,°°° er/2-2) Oa, Bn ata Pig 
L/2—-2 
t Ii 
Cope ee ee T) MET (A.7) 
i 


These calculations can be done very efficiently and clearly using a graphical 
method. 

Like in [21], one could ask for constraints on t so that the wavefunction behaves 
like the ground state of a CFT. This will be explored in further detail in a future 
work. 


B Symmetries of the Bruhat-Tits tree 


In the section 4.2 it has been discussed how to build a tensor network given a partic- 
ular graph. The global isometry of the tensor network is closely related to the graph 
symmetry. 

As already discussed in this paper recovering bulk isometry from the tensor net- 
work is an essential part establishing the bulk/boundary correspondence. In the 
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following, we will give a complete list of how the SL(2,Q,) acts on the BT tree. The 
nodes on the tree are parametrized using two vectors put into a standard form dis- 
cussed in section 5. According Iwasawa decomposition, an SL(2, Qp) group element 
G can be decomposed as 

G=NDK (B.1) 
where N is upper triangular matrix, D is diagonal scaling transformation and K is 
SL(2, Zp) transformation [80]. We note J,R € K. We explain each of them in the 

1. Dilatation 


following. 
=) 0) Cy — eR). Oe) 


which changes the accuracy by p2”. 
2. Rotation 


n= (03) (Cn + E ws 
where Jal, aie preserves the accuracy. 
3. Translation 
r=). CP a EE) ws 
which also preserves the accuracy. 
4. Reflection 
ER (6). (GF). e 


where N is the evaluation of x”) in equation (5.28). 


C Recovering SL(2, R) transformation from Coxeter lattice: 
Some examples 


In appendix B, we discussed how the SL(2,Q,) symmetry acts on the BT tree on 
which a tensor network can be defined. The BT tree was studied given the huge 
amount of symmetry it enjoys. However, ultimately the tensor networks had the 
goal of simulating conformal theories that live in the real space. In this appendix, we 
will discuss to some detail examples of networks that do preserve some subgroups of 


SL(2, R) by embedding them in regular tessellations of the AdS space. We will con- 


strast these examples with the BT tree. We will consider in particular two different 
examples where the graph is embedded in the AdS space. 
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C.1 Hexagons as a regular tessellation in the upper half plane 


As mentioned above, one way to encode part of the conformal symmetry in the 
network is to consider embedding the network in a regular tessellation of the AdS 
space. These tessellations are describable by the so called Coxeter groups, which, 
concretely in a two dimensional bulk space, is a discrete subgroup of SL(2,R). For 
example, the triangle group W(6, 2,4) is used in constructing an orbifold black hole 
in [21]. 


Figure 15. The point P is reflected through the edge QR to the new point Pı which is 
equivalent to boosting the point P to point Pı. We identify the scaling using this refection. 
Also we can identify the other two SL(2, R) generators, in fact linear combinations of them 
through the other two reflections through the other two edges of the hexagon as shown by 
the red arrows. 


Let us study this example [81] and demonstrate how various SL(2, R) generators 


can be recovered from the tessellation. To be precise, we consider embedding AdS3 
in flat space. It is thus a surface satisfying 


— X- X? + X? + XZ = —R’. (C.1) 


In the following, we will set R = 1. Also, we will take a time slice, such that X3 = 0. 
The resulting surface is a Euclidean AdS whose tessellation we will study. The 
smallest triangle (PQO) that generates the tessellation has its three angles given by 


re 
generate the full tessellation. We can choose without loss of generality where to place 


as shown in the figure. Reflection across the three sides of this triangle will 


this triangle. One convenient choice is the following: 


V3sinh(d) sinh(d) cs) R- [Ae _sinh(d) 


2 - 2 2 i 2 


cono} 
(C.2) 


P = {0,0,1}, o=] 


where d is the distance of point Q and R from the origin P, and that 


cosh(d) = v3. 


—57— 


Reflection across QR of the point P will take it to the point P}. 
This will generate a new point P 


A= {2v2,0,3} GE) 


One can easily convince that this reflection is equivalent to giving a boost (in the 
x — z plane) to the point P to Pı. The boost parameter (7) can be related to the 
angle of the triangle. For this case we have , 


n = cosh” '(3). (C.4) 


One can find successive planes such that reflection across those planes will generate 
the same amount of boost. Therefore, the above boost is in fact a global symmetry 
of the lattice. We can identify these transformations with generators of SL(2, R). In 


Poincare coordinate, 
_ dt? — d? + dz? 


ds? 2 (C.5) 
Scale transformation at the boundary is a global isometry in the bulk where 
E> Aĉ, z => àz (C.6) 


at the t = 0 slice. Using the relationship between the Poincare coordinates with the 
embedding coordinates, 


z 1 + wo 1 Z 1 — wọ 1 
), Xt = -(w+0), X2 = (1+ ar 


z2 7 5 = (w—w@), (C.7) 


where w = z + t, ù = x — t, we can identify the above boost with scaling symmetry. 
Now we have identified one such boost (B) , 


3 02,72 
B= 010 l. (C.8) 
2/20 3 


Next we define the matrix which are the generators of SL(2, R) group. 


1/01 1 1 1_1 
we ia) alaa) e(a) c9 
a (10 E -i 
They act on the embedding coordinates by 
t 
(e) M. (a=) ; (C.10) 
Xo X2+Xı 


Xı—iX> Xo 
B , so that a is obtained by comparing with the boost. In this case it is simply 


where M = . We can identify exp(aLo) with the action of 
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cosh(@) = 3. The two other planes generate other SL(2,R) transformations. We 
quote them below. 


cosh(q@) + sinh (ã) (GM + v3 vp) 
cosh(@) + sinh(@)(5Mi — ma). 


(C.11) 


In this case cosh(@) = v2, Lo = “44, L-1 = 4M; + $M, L-1 = $M3 — $M. We note 
that these type of identifications are not unique. 


C.2 Embedding of a tree in the AdS space 


Let us know look into another class of triangle group W (o0, 2, p + 1) and the tessel- 
lations generated by them. One important point about this class of groups is that it 
is shown in [82] one can embed a tree into tessellations generated by these groups. 
It is also pointed out [82] that different tessellations may have the same type of tree. 
Simplest example of this particular set is W (oo, 2,3). 


Figure 16. The figure demonstrated the W (o0, 2,3) tessellation of the upper half-plane. 
The shaded region by green strike is the smallest triangle with one vertex at œo and with 
two angles respectively 4 and 5. The point marked in blue has the coordinate z = 7 which 
is a degree two vertex. The adjacent black one has coordinate z = e'3 which is a degree 
three vertex. The embedding tree is identified by joining all the black vertices as shown in 
the diagram. Also the red arrows depict the movement of vertices (black ones) along the 
main branch under the action of T which moves z to z + 1. 


The tree generated by this group can be mapped to upper-half plane. W (co, 2,3). 
is isomorphic to PSL(2, Z). It has basically two generator S and S.T, where the 
matrices S and T are defined below by their action on the complex upper half plane, 


1 
S:z7--,T:zoz4+1. (C.12) 
z 
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Few comments are in order here. It’s obvious that we cannot recover scaling (z > az 
for a Æ 1) form this tessellation as it is just a PSL(2, Z) group as contrary to the above 
example. This has two fixed point 7 and e7 of degree two and three respectively. 

For a more general case of W (oo, 2, p + 1) the group is defined by the following 
two transformation, 


1 
a, and z — z+ 2 cos(0p). (C.13) 


where we define 0 For p an arbitrary prime, let’s repeat the exercise. The 


p= T PT 
angle between any two adjacent edges emerging from one vertex is 26,. Since the 


space is homogeneous, we can place the origin of the graph O at an arbitrary point 


in H. Let’s put it in zo = e”. It has (p+ 1) neighbors. Demanding that two of 
them sit at zo +1 fixes the positions of all of them z, = £k + i Yp with 


3 + cos 26, — 4cos(2k + 2)6, 


(3 + cos 20, — 2 cos(2k0,) — 2 cos(2k + 2)6,) 
i 29 i (C.14) 


Ty = cos Ôp 


2 sin 6, cot 
(1 + cos(k + 1)6,)(1 + cott T tan? (4**6,)) 


Yk = 


with k = 0,1,2,...,p. We can generate the whole tree by operating group action 
(C.13) on these nodes. Interestingly, these groups are known by mathematicians as 
Hecke group [83]. But still we do not have scaling. 


C.2.1 A comparison with the BT tree 


The tree structure that follows from the regular tessellations W (oo, 2, p+ 1) are also 
p+1 valent. Locally, it thus looks exactly like the p-adic BT tree. One might wonder 
therefore in what sense is the tree embedded in a regular tessellation in the upper 
half plane different from the BT tree. 

On the Bruhat-Tits tree, consider a global shift along the main branch. The 
transformation effected at the boundary is equivalent to a scaling. 


Z—> pZ. (C.15) 


This transformation surely does have a counter-part in the tree embedded in the 
upper half plane. In terms of the upper half plane complex coordinate o = x + iz, 
the transformation is basically o + ¢+1, whichis T. Recall that the boundary of the 
tree now resides on the real line of the upper half plane. This transformation therefore 
is perceived as a translation. In principle, this could be an artifact of our specific 
choice of embedding and that a suitable coordinate transformation might convert it 
into scaling at the boundary. However, one can readily show that there exists no 
non-singular coordinate transformation such that the T transformation o > ø + 1 is 
mapped to o + po. And this is only one example. Most of the symmetries discussed 
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in the previous appendix does not even have a counter part in terms of generators of 
the triangular group. Therefore, although we can geometrically embed the Bruhat- 


Tits tree H, into the hyperbolic plane such that it is the spanning tree of the regular 
tessellation given by the Coxeter group W|oo,2,p + 1], the SL(2,Q,) symmetry of 
the tree H, does not coincide with the discrete subgroup I’,. Therefore, if we are only 


intent on preserving the symmetry preserved by the given tessellation, the resulting 
wavefunction describes a very different state compared to that constructed from the 
BT tree, where the entire SL(2, Qp) is built in. 
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